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Introduction 


As the United States' space program heads toward the era of permanent occupation of 
the near Earth environment via the habitation of a permanent Space Station, the problem of 
protection of astronauts from the harmful effects of cosmic radiation assumes greater and 
greater importance. Of special concern are the effects of galactic cosmic rays in the form of 
relativistic nuclei. 1 ) When a relativistic nucleus impinges upon a spacecraft wall, it 
undergoes several reactions, an important one of which is fragmentation, whereby the 
nucleus breaks up into smaller pieces which subsequently decay. It is these secondary 
particles which contribute to the radiation environment inside a spacecraft. The problem of 
fragmentation of relativistic nuclei is the subject of a continuing research effort by the 
present author and others. 2 ) An alternative mechanism, which can produce the same sort of 
radiation environment as the fragmentation mechanism, is that of electromagnetic 
dissociation by which the electromagnetic field of one nucleus excites states in another 
nucleus which subsequently decay. The present author has led a study 3 ) of the significance 
of this mechanism to the galactic heavy ion breakup and has concluded that it is of major 
importance (i.e. large cross section) for the removal of a few nucleons when a relativistic 
nucleus impinges upon a spacecraft wall. Consequently, the mechanism of the 
electromagnetic dissociation must be included when one is trying to predict the radiation 
environment in the interior of a spacecraft when the exterior environment is composed of 
galactic heavy ions. 


Analysis of Weizsacker-Williams Theory 


The Weizsacker-Williams (WW) method of virtual quanta has been used extensively in 
the literature (see reference 4) to analyze the electromagnetic (EM) interactions of nuclei. 
The prevailing opinion seemed to be that WW theory worked quite well with a few 
exceptions. However, I believed that it was quite difficult to judge the validity of WW 
theory because authors generally used different input values of the minimum impact 
parameter, bmin, below which the reaction proceeds via the Strong Force. I decided that a 
better way to proceed would be to fit the calculated WW cross section to experiment 
treating b m in as an adjustable parameter. Thus if bmin turned out to have a ridiculous 
value, then it was absolutely clear that WW theory failed. Based on this work 5 ) it is clear 
that WW theory does not work for nucleon emission from 197 Au. Furthermore, there are 
problems with some other reactions as well. A possible explanation 5 ) for the 197 Au 
discrepancy is the use of the one-photon exchange approximation inherent in WW theory. 
Future studies of these discrepancies are strongly recommended. 

Also, in the course of this study I found that some rather serious wrong calculations 
had been presented in the literature by a group of authors. I have subsequently corrected 
the mistakes of those authors 6 > 7 ). 



Electric Onadnipole Interactions 

WW theory is also based completely on electric dipole (El) interactions 4 ), whereas it 
has been known for some time that electric quadrupole (E2) effects could perhaps 
contribute a substantial fraction to the cross section 4 ). However, calculations were not 
performed in detail because it was not known exactly how to break up photonuclear cross 
sections into their separate El and E2 components. 

I subsequently developed a calculational method which primarily uses experimental 
photonuclear cross sections as input. The El contribution is then defined as 

°E1 = a expt - <*E2 (!) 

where Og2 is calculated theoretically. The advantage of this procedure is that even if the E2 

parameters are not very well known, the total cross section will still be very accurate. This 
is because an underestimate in <% 2 will cause an overestimate in and the effects tend to 
cancel. 

It was found 9 ) that electric quadrupole effects are not significant for proton and neutron 
emission from 12 C, 16 0 or 18 0. However, E2 contributions 8 ’ 9 ) are substantial for neutron 
emission from 59 Co, 89 Y and 197 Au, generally leading to improved agreement between 
theory and experiment Notable disagreements occur for 139 La projectiles (1.26 GeV/N) 
where the theoretical cross sections are too big. Quadrupole effects improve the theoretical 
results for 16 0 projectiles at 60 and 200 GeV/N, although the theoretical cross sections are 
still too small. 

In general it was found that electric quadrupole effects are an important component in 
nucleus-nucleus collisions and that these effects can be calculated accurately 8 ’ 9 ). 


Parameterization of the Cross Sections 


A parameterization of nucleus-nucleus EM cross sections has been presented 
previously 10 * 11 ). However, since then a much more accurate theory which includes E2 
effects has been developed 8 * 9 ). In addition a much more extensive experimental data set is 
now available, including very heavy projectiles such as 139 La and energies as high as 200 
GeV/N for 16 0 projectiles. 

I therefore decided to re-examine this problem and try to come up with a more accurate 
and simpler parameterization 12 ). In this I succeeded. The new parameterization now 
includes E2 effects, and, in addition does nst require the numerical integration over energy 
that the old parameterization 10 * 11 ) required. 

Actually I have worked out three new parameterizations with accuracy varying 
inversely with simplicity 12 ). The first parameterization 12 ) (which includes E2 effects) is the 
most accurate and requires numerical integration over energy. The second 
parameterization 12 ) (which again includes E2) is based on photonuclear sum rules 4 ) and 
does m require numerical integration and is almost as accurate as the first 
parameterization. The third parameterization 12 ) (which does nol include E2) involves a 
very simple formula which can be evaluated on a pocket calculator, and it no longer 
involves evaluation of complicated Bessel functions. 

These parameterizations will be extremely useful in calculating the radiation 
environment inside a spacecraft. What needs to be done now is to implement them into the 
NASA Langly Research Center's computer transport codes. 



Other Work 


I was also involved in several other projects such as multiple nucleon knockout, two- 
neutron multiplicity, excitation energy calculations and alpha break-up. Details can be 
found in the articles in the publication list and the conference proceedings and NASA 
papers listed at the beginning of this report. 
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ABSTRACT 

The most important particle emission processes for electromagnetic excitations in 
nucleus-nucleus collisions are the ejection of single neutrons and protons and also pairs of 
neutrons and protons. Methods are presented for calculating two-neutron emission cross 
sections in photonuclear reactions. The results are in a form suitable for application to 
nucleus-nucleus reactions. 
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INTRODUCTION 

When cosmic rays in the form of heavy nuclei pass through spacecraft walls and 
astronauts' bodies, they undergo an interaction with the atomic nuclei in the spacecraft or 
astronauts. Of the four forces that we currently know about (strong, weak, electromagnetic 
and gravitational), the cosmic ray interaction occurs via both the strong and electromagnetic 
forces. The strong and electromagnetic (EM) cross sections are of comparable magnitude 
in some situations. Previous work (refs. 1-5) has concentrated on studying this 
electromagnetic aspect of nucleus-nucleus collisions. This study has so far considered only 
single-nucleon emission processes. A preliminary study has been made of multiple- 
nucleon emission (ref. 5) . However, the results of this study are unable to be utilized in 
heavy-ion transport codes because experimental photonuclear cross sections are used as 
inputs into the calculation of nucleus-nucleus EM cross sections. 

The aim of the present work is to show how to calculate these photonuclear cross 
sections for multiple-nucleon emission. Given analytic expressions for these cross 
sections, it will then be possible to add multiple-nucleon emission due to the EM effect into 
heavy-ion transport codes. Calculating multiple-nucleon emission effects for the EM 
interaction is much more difficult than the calculation for single-nucleon emission. Thus 
the present work will consider only the most important multiple-nucleon emission process- 
that of two-neutron emission . Other multiple-nucleon effects such as emission of two 
protons, a neutron and a proton, or an alpha particle are more strongly suppressed than 
two-neutron emission primarily due to the Coulomb barrier. This is especially true for 
heavy nuclei. In fact even single-proton emission is completely Coulomb-suppressed for 
heavy nuclei (ref. 1). 
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MULTIPLE NUCLEON EMISSION AND ARBITRARY 
ELECTROMAGNETIC MULTIPQLES 


Before proceeding to the study of two-neutron emission, we shall first place it in 
the broad context of arbitrary EM multipoles and multiple-nucleon emission. 

The total nucleus-nucleus EM absorption cross section cr^ is given by 

a EM = X °EM^ ^ 

x 

where the sum is over all possible species emitted in the heavy ion collision and Oem ( x ) is 
the nucleus-nucleus EM reaction cross section for producing a particular species x. The 
cross section oem ( x ) was given as equation (1) in ref. 1, but we now generalize it as (see 
also equation 2.1 of ref. 6 and equation 4.1 1 of ref. 7) 

a EM (x)=SJ <T”'’(E,x)N’“(E)dE (2) 

rt/ E 0 (x) 


where 7t can be either electric (f) or magnetic {M) and / is the order of multipolarity. Each 
term in the summation of equation (2) represents nuclear excitation by a particular EM 
multipole. Each different EM multipole n / causes a particular type of nuclear excitation. 
For example, the f\ photon field causes the nucleus to go into the giant electric dipole 
resonance (GDR) mode of oscillation parameterized as 




m 


1+ t® 2 - E GDR>® ¥Z ] 


(3) 


(see equation 8 of ref. 1 and equation 6 below) which would subsequently decay into 


e\ 


various channels x, with a probability g x , sometimes called the branching ratio. 


However, an E2 multipole (ref. 8) would cause the nucleus to go into a giant electric 

quadrupole resonance (GQR), the cross section of which would have a different 

£2 

parameterization from equation (3), and where various decay probabilities g x (e.g., 
neutron versus proton decay) mav also be different from g x l . Thus the photonuclear 



4 


reaction cross section for production of species x is some fraction of the total photonuclear 
absorption cross section 

(4) 

for a particular multipole. Equation (4) is a generalization of equation (7) of ref. 1, and 

means physically that a particular EM multipole tc /, causes a collective nuclear vibration 

n/ */ 
a abs (E) which can then decay into various channels via g x (E). Note that 

X s"' (E) =1 (5) 

X 


By combining equations (2) and (4) and assuming an energy independent branching 

it/ n/ 

ratio (replacing g x (E) with g x ), the nucleus-nucleus EM reaction cross section becomes 



<^(E)N* / (E)dE-£ 

Jl/ 


Jt/ 7t / 

g x °EM-abs 


( 6 ) 


where the nucleus-nucleus EM total absorption cross section is 

°EM - abs *j P) 

Eo(x) 

The above three equations are generalizations of the three equations on page 7 of ref. 1. 
(Note that in ref. 1 the absorption cross section was written as a EM . a bs (x) in equation 

(15) of ref. 1. The dependence on x came out because the threshold Eq(x) is the lower limit 
of integration. We prefer here simply to write a EM . abs because a true absorption cross 

section should not depend on x.) In the above equations, it is the photonuclear total 
absorption cross section (Tabs 71 / (E) which gets parameterized according to the particular 
nuclear multipole excitation. For example the nuclear fl GDR excitation is parameterized 
in equation (3). 
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If one assumes that the branching ratios are independent of multipolarity (which 
may receive some justification from the Bohr independence hypothesis (ref. 9)), then, 

i it 

replacing g x with g x , equation (6) becomes 

°EM® " g xX °EM-abs = g x °EM-abs ^ 

ill 

where the nucleus-nucleus EM total absorption cross section summed over all EM 
multipoles is 

a EM - abs = °EM - abs ^ 

it/ 

Equation (8) is interpreted physically as meaning that a nucleus-nucleus EM reaction has 
occurred exciting a nucleus into a superposition of multipolarities ^ abs ( suc ^ as 

it/ 

a linear combination of GDR and GQR). This superposition can be likened to the 
Compound Nucleus concept of Blatt and Weisskopf (ref. 9). This compound nucleus 
superposition then decays via multipole-independent branching ratios g x as in equation (8). 

Although equation (6) is probably more correct, it will be impossible to implement 
in practice due to the difficulty of calculating multipole-dependent branching ratios. 
Therefore in practical calculations, equation (8) will be used. This equation is entirely 
consistent with the concept of compound nucleus formation (ref. 9) and decay independent 
of the mode of formation (Bohr independence hypothesis). This is identical in spirit to the 
Abrasion- Ablation model (ref. 10). In fact, attempts (ref. 11) have been made to 
understand the Abrasion- Ablation model in terms of compound nucleus formation and 
decay using a T-matrix approach. Thus our basic equation (8) can be thought of as an 
Electromagnetic Abrasion- Ablation model. 

Finally, we put our previous studies (refs. 1-5) in the context of equations (1) and 
(8). First of all, equation (8) tells us that in order to calculate the nucleus-nucleus EM 
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reaction cross section, we should sum over all possible EM multipolarities (SO, SI, S2 ..., 
M0,M\,M2 ...)• The only multipolarity that has been studied so far is SI (ref. 1-5) 
leading to the Giant Electric Dipole Resonance. Clearly the effects of other multipolarities 
must be considered. Bertulani and Baur have concluded (ref. 6) that the electric quadrupole 
(ST) contribution to the total nucleus-nucleus EM cross section can be as much as 50 
percent of the SI contribution at 100 MeV/N and about 20 percent at 1 GeV/N dropping to 
about 10 percent at higher energies. Note that all the EM data (see ref. 1) are at high energy 
and so in comparing our theory (ref. 1) with experiment, the S2 contribution has not been 
large. Nevertheless we require the nucleus-nucleus (fragmentation and EM) theory to 
include all energies in the cosmic ray spectrum, and thus it is very important to consider 
other EM multipoles as well. Apart from multipoles other than SI, equation (1) tells us 
that we need also consider not only single-ucleon emission but multiple-nucleon emission 
as well, such as emission of np, 2n, 2p, a, nnp, npp and 3n. However, references 1-4 
have considered only the SI multipole and only single-nucleon emission. Thus there still 
remains much territory to explore, namely multiple-nucleon emission with SI excitation 
and then single-nucleon and multiple-nucleon emission for all other multipoles. The 
present paper is concerned only with two-neutron (2n) emission from S 1 excitation. 



7 


THEORY Of TWQ-NEUTRQN EMISSION 

Cucinotta et al. (ref. 5) have in fact studied the problem of multiple-nucleon 
emission (2n, 2p, p3n, 3n) with £ 1 excitations only. However, the photonuclear reaction 
cross section Oy(E,x) was obtained from experimental data, and so even though useful 
conclusions were drawn, the work of ref. 5 was incomplete and cannot be used in a 
transport code, which requires analytical expressions for a v (E, x). The aim of the 

present work therefore is to begin the study of multiple nucleon emission in nucleus- 
nucleus electric dipole (f\) excitation reactions using an analytic approach. The motivation 
for such an approach is to implement analytic cross section expressions into cosmic ray 
transport codes. There is no cosmic ray transport code in existence which includes 
anything other than single-nucleon emission for EM reactions. Because a full theory for 
multiple-nucleon emission is much more complicated than for single-nucleon emission, the 
present work will be limited to a study of two-neutron emission, which is the dominant 
multiple-nucleon contribution. 


TWO-NEUTRON MULTIPLICITY 

Fuller et al. (ref. 12, pp. 190; ref. 13, pp. 4; ref. 14, pp. 143-145) have defined 
the total photoneutron yield cross section as 

o(y, xn) s a(y,n) + 2c(y, 2n) + 3c(y, 3n) + .... (10) 

and the photoneutron cross section (i.e., the sum of cross sections in which at least one 
neutron is emitted) 

a(y, sn) ■ o(y, n) + a(y, np) + a(y, 2n) + a(y, na) + a(y, 3n) + .... (11) 


The neutron multiplicity is defined (ref. 15) as 

M(E) S £^ 
a(Y, sn) 


( 12 ) 


which for n and 2n emission only becomes (ref. 16) 
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M(E) 


o(Y t n) + 2o(y, 2n) 
ct(y, n) + ct(y, 2n) 


(13) 


This expression for the multiplicity has been implicitly assumed also in Thies and Spicer 
(ref. 17, equation 5). Note the reason for writing the neutron multiplicity this way: for 
c(y, 2n) equal to zero (meaning that only single neutrons (n) are emitted), the multiplicity is 
1 as expected, and for o(Y, n) being zero (meaning only two-neutron pairs (2n) are 


emitted), the multiplicity is 2 as expected. Often instead of using multiplicity, one works 


with the quantity 

q(Y> 2n) 

a(Y, n) + o(Y, 2n) 


= M(E) -1 


(14) 


Blatt and Weisskopf (ref. 9, chapter VIE, section 6B) have worked out the general theory 
of multiple-nucleon emission, which they call secondary nuclear reactions. They write the 
cross section for production of particle b (see ref. 9, equation on pp. 375), which we 
specialize to photoproduction as 

o(y, sb) = a(Y, b) + ^ a(y, be) (15) 

c 


This is made up of a primary, single-step cross section o(y, b) in which particle b is emitted 

directly from decay of the compound nucleus and a secondary, multiple-step cross section 
o(y, be) is which either particle b or c is emitted from the compound nucleus which decays 
to a lower excited state which again decays via emission of particle c or b. Blatt and 
Weisskopf (ref. 9, pp. 376-379) then go on to calculate a(y, b) and a(y, be) in terms of 
a(y, sb); however, one generally needs to work the integrals out numerically. We shall 
eventually do this when working out the general expression for multiple-nucleon emission. 
However, for the case of two-neutron emission, certain simplifying assumptions 
concerning the integrals can be made so that analytic expressions can be obtained. These 


are (ref. 9, eqn. 6.14, pp. 377) 

o(y, n) = a(Y, sn) (1 + Esec/©) exp(-£sec/©) 


(16) 
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and (ref. 9, eqn 6.18, pp.379) 

a(Y, 2n) * a(y, sn) [1 - (1 + esec/Q^xpC-Csec/©)] (17) 

where 0 is the nuclear temperature and 

£sec = Ey - E 0 (Y, 2n) (18) 

The above equations are also discussed in references 15, 16, 17, 18. Equation (17) is 
derived on the assumption that the photoneutrons are produced via statistical decay of the 
compound nucleus (ref. 9). However this is not always necessarily true. It can sometimes 
happen that the incident photon, rather than exciting a compound nucleus which 
subsequently decays via neutron emission, will rather knock a neutron out directly. (Note 
that this direct emission violates the abrasion-ablation concept.) We therefore introduce a 
fraction of direct emission f d , which typically has values (ref. 15, and G. O'Keefe and R. 
Rassool, private communication) in the range 0.1 - 0.2 although it can be as large as 0.4 
(ref. 17). This direct fraction is incorporated into our theory (ref. 15, 17) by generalizing 
equation (17) as 

c(y, 2n) * a(y, sn) (1 - f d ) [1 - (1 + e se c/©) exp (-£<«*/©)] (19) 

This equation is our main result for 2n emission. The quantities f d , £sec> and © 316 easily 
calculated so that the only input required is a(y, sn). Blatt and Weisskopf (ref. 9, pp. 379) 
claim that a(y, sn) can usually be closely approximated by the compound nucleus formation 
cross section, i.e., the total absorption cross section. (This is done in reference 18.) 
However, this is only true when compound nucleus decay proceeds predominantly by 
neutron emission which is the case for heavy nuclei. For light nuclei, where proton 
emission is just as important (ref. 1), the above approximation is not valid. Thus we now 
neglect all contributions to o(y, sn) other than a(y, n) and a(y, 2n) so that 

a(y, n) + q(y, 2n) ' (1 ‘ f d> + 

which is the result written in reference 16 and 18 for f d = 0. Defining 
Xffac = (1 - f d ) [1 - (1 + £ S ec/@)exp(-e sec /0)] 


( 20 ) 

( 21 ) 
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we rearrange (20) to obtain 
Xfrac 

° ( * 2n) = TTxf^° ( * n) ' (22) 

Although equation (19) is our fundamental equation for 2n emission, we shall in fact use 
equation (22) to test the general theory. This is because several authors have reported both 
o(y, n) and a(y, 2n). We shall use experimental input for o(y, n) and then calculate 
a(y, 2n) using (22) and see how well our calculations agree with experimental 
measurements. If we obtain good results, we can be confident that the theory outlined 
above can be used to calculate 2n emission in EM nucleus-nucleus collisions. 

Finally, we note that the above equations can be used for an approximate calculation 
of the neutron multiplicity given in equation (13) as 

M(E) = 1 + (1 - fd)[l - (1 + e se c/©)exp(-e S ec/©)] (23) 

which is the result reported in references 15 and 17. (Note however the error in equation 
(3.12) of reference 15, which is that cr(y, np) should not appear in the numerator). 
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NUCLEAR TEMPERATURE 

In the above expressions, the nuclear temperature 0 plays a central role. For a 
Maxwell-Boltzmann (classical) gas (ref. 19, pp. 1 17), the relation between energy and 
temperature is linear 

E = |n k© (24) 


whereas for a Fermi gas, the relation is quadratic 

E = (const.) k 2 0 2 (25) 

The theory developed by Blatt and Weisskopf (ref. 9, pp. 372) uses a Fermi gas nuclear 
temperature given by 



where the constant A depends on the total nucleon number A. We shall not use the Blatt 
and Weisskopf theory for nuclear temperature, but rather use a more sophisticated version 
developed by Bohr and Mottelson where (ref. 20, eqn. 2-50, pp. 154) 

e'=-L^ 


p BE 

with the Fermi gas energy level density (ref. 20, eqn. 2-47, pp. 153) 

6 1/4 g o \ 1/2 

p(N, Z, E) = — exp^ g e) ] N » Z 

(g 0 E ) 5/4 


(27) 


(28) 


being of similar form to that of Blatt and Weisskopf (ref. 9, eqn. 2-50, pp. 154). The 
symbol E is the excitation energy, and g 0 is the "one-particle level density at the Fermi 
energy, representing the sum of neutron and proton level densities" (ref. 20, pp. 153). 
Thus the nuclear temperature becomes (ref. 20, eqn. 2-50, pp. 154) 



12 


which is the same form as used in reference 15. It is very important to realize that the 
energy E in the above expressions is, in fact 


E — Ey - Eo(Y, n ) 


(30) 


given by the single-neutron emission threshold Eq(y, n). This is emphasized in references 

18 (pp.304) and by G. O'Keefe and R. Rassool (private communication). The one-particle 

level density go is related to the level density parameter a via (ref. 20, eqn. 2-123, pp. 187) 
2 



o’ 


(31) 


which is plotted in Figures 2-12 of reference 20 (pp. 187). It is seen that a very much 
depends on shell structure, and therefore whenever possible, a should be obtained directly 
from Figures 2-12 (ref. 20). However, a rough approximation is 

a-ylMeV' 1 ) (32) 


but note that this fails badly for A between 190 and 210, (ref. 20). 
Rewriting the nuclear temperature as 



(33) 


with the first term corresponding to the form used by Blatt and Weisskopf (ref. 9), we 
approximate a to obtain 
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TEST OF TWO-NEUTRON MULTIPLICITY THEORY 

A computer code was written to implement the theory described above. Necessary 
inputs are the energy level density parameter a which was taken directly from Figs. 2-12 
(ref. 20, pp. 187) rather than from equation (32). The two-neutron cross section was then 
calculated from the single-neutron cross section with a best-fit value of fd, the fraction of 
direct emission. The results are listed in Tables 1-8. As can be seen, the agreement 
between theory and experiment is good, except for a few notable exceptions. It is 
disturbing, however, that the direct emission fractions are so high. 


% 
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DISCUSSION 

We have demonstrated that the theory of two-neutron emission described herein is 
able to predict most two-neutron multiplicities with a high degree of accuracy. One 
therefore wishes to implement this theory into the calculation of nucleus-nucleus cross 
sections and eventually into a transport code. A parameterization of the energy level 
density parameter a can be very easily worked out (ref. 20). Thus the only unknowns are 
fd and o(y, n). 

The values of fd were obtained from an overall best fit to c(y, 2n) data. (In 
practice, a single o(y, 2n) value at a single energy determines fd.) Clearly some new work 
remains to be done here. One can either work out a parameterization by determining fd 
from best fits to a whole range of data, or better still, one should determine fd theoretically. 
Our initial effort will undoubtedly be a parameterization, so that the theory can be used 
within transport codes. 

The most serious unknown is c(y, n) which is used as input for the o(y, 2n) 
calculation at each energy. One might think that the liquid drop model calculations for 
a(y, n), described extensively in reference 1, would be adequate for this purpose, 
especially as these same model calculations gave such good results for nucleus-nucleus 
cross sections. Unfortunately, this is not the case. Even though the liquid drop model 
calculations of a(Y, n) are quite accurate at most photon energies, they are unfortunately not 
very accurate for high photon energies. The reason that this did not matter for the 
calculations of reference 1 was that the virtual photon spectrum N(E) predominated at low 
photon energies, and so energy integrated nucleus-nucleus cross sections were quite 
accurate. However, obviously (y, 2n) processes only occur in the high energy (y, n) 
region. The o(y, 2n) calculations require a(y, n) as input, and thus we now need very 
accurate calculations for a(Y, n) at high energy. (Calculations described in the present 
paper used experimental numbers.) 
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Thus, the most pressing need for inclusion of the effects of multiple-nucleon 
emission in heavy ion transport codes is the accurate calculation of a(y, n) at high energy. 
The first such step will require the study of deformation splitting (ref. 14, 17, 18, 24, 25) 
and isospin splitting of the giant dipole resonance (ref. 24, 25). 

Finally, even though only two-neutron emission has been considered herein, it is 
now clear how to proceed with the general multiple-nucleon emission problem. One can 
simply numerically integrate the expressions of Blatt and Weisskopf (ref. 9) to obtain any 
multiple-nucleon final state. (These expressions simplified in the two-neutron case 
considered herein, so that numerical integration was not necessary.) This method has 
previously been applied to 32 S(y, d) by Norbury (ref. 26). 

As an alternative to the whole approach taken in the present work (and even for the 
single-nucleon work of ref. 1), it may be feasible to instead do calculations using some of 
the modem nuclear evaporation codes (e.g., EVAP-4 and EVA-3). The advantage of these 
codes is that they fully incorporate both statistical and direct emission. The excitation 
energy in these codes is simply the photon energy considered above. 

CONCLUSIONS 

Electromagnetic excitations is nucleus-nucleus collisions can often involve the 
ejection of pairs of nucleons. Methods have been presented for calculating two-neutron 
emission cross sections in photonuclear reactions. These cross sections are now in a form 
suitable for use in nucleus-nucleus reaction theory. 
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Table 1 12 C photoneutron cross sections. Data are from ref. 21 


Ey, MeV 

c(y, n), mb 

a(y, 2n), 
Experimental 

mb 

Theoretical 
fd = 0.9 
a = 2 MeV -1 

35 

2 

0.004 

0.043 

40 

0.8 

0.012 

0.051 

50 

0.2 

0.015 

0.020 

60 

0.1 

0.014 

0.011 

80 

0.1 

0.095 

0.011 
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Table 2 18 0 photoneutron cross sections. Data are from ref. 27 


Ey, MeV c(y, n), mb a(y, 2n), 


Experimental 


mb 


Theoretical 
f d = 0.5 
a = 3 MeV' 1 


14 

18 

22 

26 

30 

34 

38 

42 


4 

2.4 

4 
7 

5 
2 
2 
1 


1 

3 

3 

3 

3 

3 

2 

3 


0.5 

1.3 

3.1 
6 

5 

1.9 

1.9 

1.2 
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Table 3 40 Ar photoneutron cross sections. Data are from ref. 22 


Ey, MeV a(y, n), mb o(y, 2n), 


Experimental 


mb 


Theoretical 
f d = 0.05 
a = 6MeV-l 


18 

20 

21 

22 

23 


16 

7 

5 

3 

2 


5 

14 

13 

11 

10 


6 

12 

14 

12 


11 
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Table 4 59 Co photoneutron cross sections. Data are from ref. 28 


Ey, MeV o(y, n), mb o(y, 2n), mb 


Experimental Theoretical 

f d = 0.6 
a = 8 MeV- 1 


18 

65 

0.3 

27 

20 

47 

9 

4 

22 

21 

17 

7 

24 

17 

12 

9 

26 

14 

11 

8 

28 

19 

9 

12 

30 

20 

5 

13 

32 

9 

6 

6 

34 

1 

6 

1 


3 3 


36 


4 
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Table 5 63 Cu photoneutron cross sections. Data are from ref. 23 


Ey, MeV a(y, n), mb o(y, 2n), mb 


Experimental Theoretical 

fd = 0.3 
a = 8 MeV' 1 


20 

44 

2 

0.6 

21 

26 

4 

5 

22 

20 

8 

11 

23 

13 

12 

12 

24 

11 

13 

13 

25 

10 

13 

15 
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Table 6 89 Y photoneutron cross sections. Data are from ref. 29 


Ey, MeV 

o(y, n), mb 

o(Y, 2n), 
Experimental 

mb 

Theoretical 
fd = 0.5 
a = 12 MeV -1 

20 

70 

0 

36 

22 

35 

7 

6 

24 

10 

10 

6 

26 

2 

16 

2 

28 

0.5 

16 

0.5 
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Table 7 197 Au photoneutron cross sections. Data are from ref. 16 but have been 
multiplied by 0.93 following the suggestion of ref. 30 


MeV 

a(y, n), mb 

a(y, 2n), 

mb 



Experimental 

Theoretical 
fd = 0.4 
a = 24 MeV' 

15 

391 

9 

20 

16 

214 

73 

121 

17 

93 

98 

98 

18 

56 

93 

74 

19 

37 

65 

53 

20 

27 

58 

40 

21 

17 

42 

25 

22 

11 

47 

16 

23 

7 

40 

11 


24 2 


37 


3 
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Table 8 208 Pb photoneutron cross sections. Data are from ref. 16 


Ey, MeV o(y, n), mb a(y, 2n), mb 


Experimental Theoretical 

fd = 0.2 
a = 8 MeV -1 

16 180 80 104 

18 30 80 51 


20 


20 


50 


55 
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LIST or SYMBOLS 

A nuclear temperature constant, MeV' 1 

a energy level density parameter, MeV' 1 

E energy, MeV 

EM abbreviation for "electromagnetic" 

Egdr giant dipole resonance central energy, MeV 

Eq(x) photonuclear reaction threshold for a particular species x, MeV 

Ey photon energy, MeV 

£ / electric multipole of multipolarity / 

fd fraction of direct, non-statistical emission 

go one particle level density at the Fermi energy, MeV' 1 

n t 

g x energy independent branching ratio for species x due to an/ nuclear excitation 

n/ 

g x (E) energy dependent branching ratio for species x due to a 7t /nuclear excitation 

k Boltzmann constant 1.3805 x 10' 23 J/K 

1 multipolarity 

M / magnetic multipole of multipolarity / 

mb millibam 

N 0 number of particles 

N n ^(E) virtual photon number spectrum of a particular multipole, MeV' 1 
n neutron 

p proton 

x particular particle species emitted in a heavy ion collision 

a alpha particle 

T giant dipole resonance width, MeV 

£sec maximum energy available for emission of secondary particles, MeV 
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n f particular electromagnetic multipole (e.g., £\, tfO) 
p energy level density, MeV -1 

a cross section, mb 

a(y, sb) sum of cross sections in which at least one particle b is emitted 

o(y, sn) photoneutron cross section, mb 

a(y, xn) total photoneutron yield cross section, mb 

71 1 

a abs(E) photonuclear total absorption cross section, mb 

a EM total nucleus-nucleus EM absorption cross section, mb 

<t E m( x ) nucleus-nucleus EM reaction cross section for production of a particular species 
x, mb 

CT EM-abs nucleus-nucleus EM total absorption cross section summed over all EM 
multipoles, mb 

K t 

a EM - abs nucleus-nucleus EM total absorption cross section for a particular multipole, mb 
o m giant dipole resonance cross section peak value, mb 

7t t 

°v (E, x) photonuclear reaction cross section for production of a particular species x due 
to nuclear excitation from a particular multipole n /, mb 
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Abstract 

The invariance of Classical Electromagnetism under CPT is studied by considering 
the motion of a charged particle in electric and magnetic fields. Upon applying CPT 
transformations to various physical quantities and noting that the motion still behaves 
physically demonstrates invariance. 


1. Introduction 


In teaching a recent course in electromagnetism I discussed the transformation 
properties of the electric E , and magnetic fields B under parity (P) time reversal (T) and 
charge conjugation (C) transformations. I was amazed to find that a unified discussion of 
CPT transformations in classical electromagnetism does not exist in any of the standard 
text books. The best I could find was a discussion of time-reversal and axial and polar 
vectors in (Jackson 1975). The aim of the present paper is purely pedagogical. I wish to 
discuss CPT transformations, not from a mathematical point of view (Rosen 1973) but 
from a completely physical point of view by considering how the transformations affect the 
motion of a charged particle in an electromagnetic field. This approach was very successful 
with the students; they came away with an excellent physical understanding of these 
transformations. 

Our procedure will be as follows: we shall look at the motion of a charged particle 
— ^ ^ 

in either E or B fields and then apply one of the C, P or T transformations to the whole 
situation. Each quantity, such as force, charge, etc. will transform in their own particular 
way. Because Maxwell's equations are invariant however we know that our transformed 
situation shouldn't be any different. In particular, bv looking simply at trajectory of the 
particle, we shouldn't be able to tell that a transformation has taken place . Examining these 
"pictorial" transformations will give us a very clear physical understanding of the 
corresponding mathematical transformations. 

Before proceeding we note the transformation properties of some of the quantities 
we shall need. Details can be found in (Jackson 1975, Perkins 1987, Rosen 1973). 

Electric charge q is a scalar under P and T, but a pseudoscalar under C. (Pq = q, Tq = q, 

— > — > 

Cq = -q) Velocity v changes sign (vector) under T because of its form as d x /dt. 

Otherwise it is a vector under P and a pseudovector under C. (Pv = -v,Tv = -v , Cv = v ) 

Force F or acceleration is a pseudo v ector under T because of the form d 2 x /dt 2 . Otherwise 

— > — > — > — > — > — > — > 
it transforms the same as v . (PF = -F , TF = F , CF = F ). 


2. Electric Field 

— > — > 

The transformation properties of the electric field are PE = -E (vector), 

— > — > — ► — > 

TE = E (pseudo- vector), CE = -E (vector). Let us study these in turn by considering the 
motion of a positively charged particle in an electric field, as shown in Fig. 1. 


2.1 Parity 

Parity transformations are most easily visualized by reflection in a mirror. Imagine 

that the charge in Fig. 1 is moving to the right at right angles toward a mirror. Let us now 

reflect this whole situation in the minor. Under reflection charge q is a scalar and electric 

field E , force F and velocity v are vectors and so reverse sign as shown in Fig. 2. Thus 

— ^ > 

the relative orientation of F and E are still the same. Thus we still have acceleration 
towards the mirror: we cannot tell that we are now in a reflected world. 
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Charge Conjugation 


— ► -> — > 

Under charge conjugation F and v remains the same, but E and q change sign, as 
shown in Figure 3. Thus we still have acceleration of the charge to the right, so we can't 
tell that a C transformation has been made. 



2J Time-Reversal 



F remains the same also. This time F and v become o pposite to each other. As depicted 

in Fig. 4, the motion is like a movie of Fig. 1 being run backwards. Thus the particle starts 

— ^ 

at the far right hand side with velocity -v but slows down because the force is to the right, 
as it travels into the past. The motion still makes physical sense and we cannot distinguish 
travel into the future from travel into the past. (The position of the particle in Fig. 4 can be 
superposed directly onto Fig. 1) 


3. Magnetic Field 


— > — > — > — > — > 

The magnetic field B transforms as PB = B (pseudovector), TB = -B (vector), 

— > — > 

CB = -B (vector). We consider a positive charge moving counterclockwise in a magnetic 
field as shown in Fig. 5. 


2JL Parity 

— > — > 

We now imagine a mirror parallel to the plane of the paper. Because F and v are in 

the same plane they will not be affected this time by a parity transformation. Any vector 

— ) 

perpendicular to the plane of the paper will change sign. However, B (which is 
perpendicular) is a pseudo- vector under P and will not change sign. Thus the mirror- 
reflected version of Fig. 5 is identical to that figure. Obviously then the motion is invariant 
under reflection. 


3.2 Char ge Conjugation 

Figure 6 represents a charge conjugated version of Fig. 5. As can be seen the 
motion is still counterclockwise, so we couldn't tell that a transformation has been made. 

2i 3 Time-Reversal 

The time-reversed situation of Fig. 5 is shown in Fig. 7. The motion here is 
identical to what we would see if we took a movie of the motion of Fig. 5 and then ran it 
backwards. 



Summary 


Note the importance of what has been done in the 6 transformations presented 
above. We have transformed each quantity separately, thus making it quite possible to 
obtain a resultant unphvsical trajectory, ie. a trajectory violating Maxwell's equations and 
the Lorentz force law. (An example of an unphysical trajectory is a positive charged 
particle moving clockwise in a magnetic field directed into the page.) However, none of 
our transformed trajectories have been unphysical indicating invariance of Classical 
Electromagnetism and thus preventing us from assigning an absolute value to C, P, or T. 
That is, by examining the physical motion we cannot tell if a C, P or T transformation has 
been made. It is hoped that the ideas presented in this paper will provide a useful physical 
supplement to the study of classical electromagnetism. 
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Figure Captions 


Fig. 1 . Motion of a Positive Charge moving towards a mirror with an initial velocity v in a 
Uniform Electric Field. TTie position of the charge at various times is indicated with a 

symbol ®, at time ti. 


Fig, 2 Motion of Fig. 1 reflected in a mirror 

— ► — > — > — ► — > — > 

( PE = -E, Pq = q, Pv = -v, PF = -F) 

Fig. 3 Charge - conjugated version of Fig. 1 . The positive charge becomes a negative 
charge indicated with a symbol . 

— ^ ^ — > — ¥ — > — > 

(CE = -E, Cq = -q, Cv = v , CF = F 

Fig. 4 The charge of Fig. 1 travelling into the past. 

— ► — ► — > — > — > — > 

(TE = +E, Tq = q, Tv = -v, TF = F) 


Fig. 5 Motion of a Positive Charge with a velocity v moving counterclockwise in a 
Uniform Magnetic Field directed into the page. 

— ^ ^ 

Fig, 6 Charge-conjugated version of Fig. 5 ( B points out of the page) (CB = -B) 

Fig. 7 The charge of Fig. 5 travelling into the past. 


(TB =-B) 
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ABSTRACT 

An optical model description of momentum transfer in relativistic heavy ion 
collisions, based upon composite particle multiple scattering theory, is presented. 

The imaginary component of the complex momentum transfer, which comes from 
the absorptive part of the optical potential, is identified as the longitudinal 
momentum downshift of the projectile. Predictions of fragment momentum 
distribution observables are made and compared with experimental data. Use of the 
model as a tool for estimating collision impact parameters is discussed. 

1. INTRODUCTION 

Since the pioneering experiments on relativistic heavy ion fragmentation 
using carbon and oxygen beams , 1 ' 2 attention has been directed toward 
understanding the underlying mechanisms of fragmentation processes. Over the 
past two decades, a substantial body of literature has resulted from studies of these 
phenomena, and several excellent reviews have been written . 3 " 6 Perhaps the most 
significant findings of the early experiments were the observations that the 
fragment momentum distributions were Gaussian in the projectile rest frame, and 
that the isotopic production cross sections factored into a product of target and beam- 
fragment terms. Initial attempts to explain these phenomena utilized a statistical 
model to describe the reactions . 7 ' 9 This later evolved into a two-step model called 
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abrasion-ablation 10 where the abrasion stage can be formulated using geometric 10 ' 11 
or quantum mechanical arguments. 12 ' 13 In the present work, we use the impulsive 
excitation energy ideas of Fricke, 14 within the context of composite particle multiple 
scattering theory, to derive a method for predicting momentum transfers occurring 
in relativistic heavy ion collisions. This momentum transfer is a function of impact 
parameter. A new feature of this work is that the momentum transfer is a complex 
quantity. The real component is the usual transverse momentum transfer resulting 
from elastic scattering. The imaginary component is explicitly shown to be the 
longitudinal momentum transfer, or downshift, arising from the absorptive part of 
the complex optical potential. Using this formalism, projectile nucleus fragment 
momentum "downshifts" resulting from the dynamics of the nuclear collision can 
be calculated and compared with laboratory beam measurements. In addition, 
modifications to the widths of the momentum distributions can be estimated using 
the formalism. 

The outline of the paper is as follows. In section 2 the dynamical momentum 
transfer expression is derived, and representative calculations of momentum 
transfer as a function of impact parameter are presented. In section 3 the 
connections between collisional momentum transfer and fragment momentum 
downshifts/ widths are made. A method of choosing appropriate impact parameters 
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for each fragmentation channel is then described. Next, calculations of momentum 
downshifts for fragments produced by oxygen nuclei colliding with various targets 
are made and compared with experimental data. 2 We also compute widths of 
momentum distributions for 139 La fragments and compare with recent experimental 
measurements. 16 In section 4 we propose a method for using the momentum 
transfer model to estimate collision impact parameters. Finally, in section 5 we 
conclude by summarizing the current status of model development and discuss 
future directions for research. 

2. METHOD OF CALCULATION 

In reference 17, a coupled-channels Schrodinger equation for composite 
particle scattering, which relates the entrance channel to all of the excited states of 
the target and projectile, was derived by assuming large incident projectile kinetic 
energies and closure of the accessible eigenstates. The equation is written as 

(v 2 + k 2 ) (x) = 2m A p A T (A p + A T ) * l £ V np> n y (x) y nV (x) 

t ' 

nH 

where the subscripts n and p (primed and unprimed) label the projectile and target 
eigenstates; m is the nucleon mass; A p and A T are the mass numbers of the 
projectile and target; £ is the projectile momentum relative to the center of mass; 
and x is the projectile position vector relative to the target. In terms of the nucleon- 
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nucleon scattering t-matrix t a j, and the internal state vectors of the projectile gP ^ ^ 
and target it was also demonstrated that the potential matrix is expressible as 

Vnjx, n y (x) = (gS gj V opt t Pi i T , x 


P T 
gn'gp' 


( 2 ) 


where 


^opt j^P, ^T, x J — S ^otj 

<*j 


( 3 ) 


This same formalism can be used to investigate heavy ion collision momentum 
transfers. Within the context of eikonal scattering theory, the solution to the 
Schrodinger equation 

H1'(x,£p,4 T ) = E4'(x,£p,£ r ) 
is 



1 

t 

r 

N 

( x > ^P, ^t) = ( 2lt )~ 2 eX P 

* y J V opt ( x » £p. £tJ dz 


( 5 ) 


SK «p> gj (St) <= ik 1 


where v is the velocity. The total momentum of the projectile is then given by the 
matrix element involving the sum of the projectile single-nucleon momentum 
operators as 



'F 


a p - 

E v P , 0 


a= 1 



( 6 ) 


where the subscript P on the gradient operator denotes that the gradient is to be 
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taken with respect to the projectile internal coordinates. Equation (6) actually 
denotes a potential matrix p n ^ n y in analogy with (2). Therefore, substituting (5) 
into (6) yields 


IV, n v 


= ( 6n (^p) gji (£t) 



f A P - ^ 


e -iS 

-i 1 v P , 0 

1 ) 

e iS 


(7) 


where 

S * J:J V opt (x, ^ P , ^ x )dz'. 

— oo 


( 8 ) 


Equation (7) can be further expressed as 


Pnji, n'n' ~ Po + ( 8n 


/ Ap 

-S V P.aS 


A A, 


(9) 


where the momentum before the collision is 


Po=( g£gl 


A P — 

-i I Vp.a 

a = 1 




( 10 ) 


The total momentum transfer to the projectile is then given by 


Q ~ Pn|i, n'n' Pq ” ( gn gj 


_S Vp, a S 


a 


P T 

gn'Sj' 


(ID 


For scattering near the forward directions, the couplings between various excited 
states is small and the off-diagonal elements in Eq. (11) can be neglected; hence, the 
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momentum transfer can be approximated by 


Q-Pei-Po = 


gogl 


v P>a s 



( 12 ) 


In terms of projectile and target number densities, and the constituent-averaged two- 
nucleon transition amplitude, 18 T, Eq. (12) becomes 

Q (b) = -ApA T | d 3 £p pp (^ P ) | d 3 py (£ t ) 


v P )2;(b + ?+4p-i T )^ 


(13) 


where the integration limit in the longitudinal direction has been extended to 
infinity. The momentum transfer in (13) is therefore only a function of the impact 
parameter of the collision. The projectile and target number densities (p p and p T ) 
are normalized to unity as 


j p (x) d 3 x = 1. ^ 

The constituent-averaged two-nucleon transition amplitude is obtained from the 

first-order t-matrix used in our previous studies 13 of nucleus-nucleus collisions as 


t (e, x) = - ( e / m f 2 a (e) [a (e) + i ] [2xcB(e)] ** exp [ - x 2 /2B(e)] 


(15) 


where e is the two-nucleon kinetic energy in their center-of-mass frame, o(e) is the 
nucleon-nucleon total cross section, a(e) is the ratio of the real-to-imaginary part of 
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the forward scattering amplitude, and B(e) is the nucleon-nucleon slope parameter. 
Values for these parameters taken from various compilations are listed in 
reference 18. 

The dynamical momentum transfer to the projectile, given by Eq. (13), results 
from interactions with the target. A new feature, unique to this work, is that it is a 
complex quantity. The real part of the momentum transfer, which comes from the 
real part of the complex optical potential, is the contribution arising from elastic 
scattering. It is purely transverse. The imaginary component, which comes from 
the absorptive part of the complex optical potential, is the longitudinal kinetic 
momentum downshift. To demonstrate this last assertion, we rewrite Eq. (13) 
symbolically as 

Q = (Q R + iQi)b (16) 

where i is V- 1 and £ is the unit vector transverse to the beam direction. If z denotes 
the unit vector in the incident beam direction, then from elementary complex 
analysis 19 we know that 

ib = -z (17) 

since i is an operator which rotates a unit vector counterclockwise through rc /2 


radians. Therefore, the momentum transfer is 
Q = Q R b-Q,z 


( 18 ) 
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which we relabel for clarity as 
Q = Qj.b-Q,,z. 

Note the similarity of this argument to that of complex indices of refraction in 
electromagnetic wave propagation. From Eq. (13), the transverse component is 

= ~ ApAi> | d 3 £p pp (£ P ) | d 3 ^7 p T (£j) V P J^Re t (b + z' + Ijp - £ T ) — 


( 19 ) 


( 20 ) 


and the longitudinal component is 


Qj j = - ApA T j d 3 ^ P p P (^p)j d 3 Px (^ T ) V p Im t (b + z' + - £ T ) 


( 21 ) 


Calculated momentum transfers obtained using equations (20) and (21) are 
displayed in Figure 1 for 16 0 at 2.1 AGeV colliding with a beryllium target. These 
calculations utilize the harmonic well nuclear densities from our previous 
work . 13 ' 18 From the figure, two features are readily apparent. First, the longitudinal 
momentum transfer is larger than the transverse indicating the primarily 
absorptive nature of the nuclear collision at this energy. Second, the predicted 
momentum transfers decrease rapidly with increasing impact parameter. This will 
be a subject of further discussion in subsequent sections of this paper. 
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3. RESULTS 

The collisional momentum transfers computed using the model described in 
the previous section can be related to experimentally-measured, heavy-ion fragment 
momentum downshifts/widths through considerations of energy and momentum 
conservation. As has been pointed out elsewhere, 8 ' 20 a momentum transfer in any 
direction Qj modifies the width hj of the momentum distribution in that direction 

by 



( 22 ) 


and the mean by 

p/-P j+ £Q, (23) 

From the latter, the longitudinal momentum downshift is given by 

4p ii = p lr p ii=f < 2ii <24) 

where Q ( ( is the magnitude of the longitudinal momentum transfer [obtained from 
eq. (21)], F is the fragment mass number, and A is the initial mass number of the 
fragmenting nucleus. Recalling that Q ( ( is a function of impact parameter, an 
appropriate method for choosing it for each fragmentation channel is necessary. 
Recently a semiempirical abrasion-ablation fragmentation model, NUCFRAG, was 
proposed. 21 Although it assumes simple uniform density distributions for the 
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colliding ions, and a zero-range (delta function) interaction, it does include frictional- 
spectator-interactions (FSI) and agrees with experimental cross section data to the 
extent that they agree among themselves. Also, and most importantly for this work, 
it is easily modified to yield the impact parameters for each fragmentation channel. 
Hence, the procedure for evaluating equations (22) and (24) is to extract impact 
parameters from NUCFRAG for each nucleon removal corresponding to AA = 1, 2, 

3, These impact parameters are then inserted into Eqs. (20) and (21) to obtain 

the corresponding momentum transfers for use in evaluating Eqs. (22) and (24). 
Because NUCFRAG uses uniform densities, uniform densities are also used in 
evaluating (20) and (21). In addition, the zero-range interaction in NUCFRAG is 
simulated for numerical integration purposes in (20) and (21) through the use of a 
very narrow Gaussian form for the t-matrix given by eq. (15). This narrow Gaussian 
is the same width for all collision pairs and therefore is not an arbitrarily adjusted 
parameter. 

Representative calculations for momentum downshifts as a function of 
fragment mass number are displayed in Figure 2 for 16 0 projectiles at 2.1 AGeV 
colliding with targets of Be, C, Al, Cu, Ag, and Pb. These momentum downshifts are 
target-averaged using simple arithmetic averaging. For comparison, the target- 
averaged experimental data from reference 2 are also also displayed. For display and 



12 


comparison purposes, the latter are also averaged over all isotopes contributing to 
each fragment mass number using 


5>iK 




expt 


1 °, 


( 25 ) 


where Oj is the experimental production cross section for the ith fragment isotope. 
Comparing the theoretical estimates to the experimental data, reasonable agreement 
is obtained considering the simplified form of the nuclear fragmentation model 
used in the calculations and the overall sensitivity of the calculated momentum 
transfer to the choice of impact parameter. Improved agreement is expected if 
impact parameters from a fragmentation model using realistic nuclear densities and 
interactions were available. This is especially true for collisions involving lighter 
ions, such as carbon, oxygen, and beryllium, which are poorly represented by simple 


uniform nuclear distributions. 

Figure 3 displays transverse momentum widths as a function of fragment 
mass number for 1.2 AGeV 139 La fragmenting in carbon targets. The experimental 
data are taken from reference 16. Again, impact parameters from NUCFRAG are 
used as inputs into the momentum transfer expressions [Eqs. (20) and (21)]. For 
consistency with the use of these impact parameters, a narrow Gaussian t-matrix 
and uniform nuclear densities were again utilized in the momentum transfer 
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calculations. From Figure 3, it is clear that the agreement is much better than in 
Figure 2 and probably reflects the fact that a uniform nuclear density distribution is a 
more reasonable aproximation for a heavy nucleus-like lanthanum than for light 
nuclei such as oxygen. 

4. ESTIMATING COLLISION IMPACT PARAMETERS 

Thus far in this work, we have used collision impact parameters as inputs 
into a momentum transfer computational model which in turn has yielded 
estimates of heavy ion fragment momentum downshifts /widths for comparison 
with experimental data. However, this procedure can be reversed and the model 
used to estimate collision impact parameters from measured momentum 
downshifts for relativistic collisions. Let F be the fragment mass number with 
measured longitudinal momentum downshift AP , ( produced in a relativistic 
collision between a projectile nucleus (mass number A) and some target. Then, 
from eq. (24), the longitudinal momentum transfer to the projectile from the target 
is 

Q| = A.AP||. (26) 

The collision impact parameter can then be estimated from eq. (21) by computing 
Qi | as a function of impact parameter (e.g., in Figure 1) and using Q, , from eq. (26) 
as the entry. To illustrate, consider a collision involving 2.1 AGeV oxygen colliding 
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with a beryllium target. The calculated momentum transfer using realistic nuclear 
densities are displayed in Figure 1. If the measured (hypothetical) momentum 
downshift for the 14 N fragment is 35 ± 7 MeV/c r then eq. (26) yields a longitudinal 
momentum downshift of 40 ± 8 MeV/c. From Figure 1, the corresponding range of 
impact parameters is 6.1 - 6.4 fm. A similar procedure incorporating measured 
momentum distribution widths and Eqs. (22) and (20) or (21) could also be used to 
estimate collision impact parameters. 

5. CONCLUDING REMARKS 

Beginning with composite particle multiple scattering theory, an optical 
model description of collision momentum transfer in relativistic heavy ion 
collisions was derived. General expressions for transverse and longitudinal 
momentum transfer, which utilize a finite-range two-nucleon interaction and 
relativistic nuclear densities, were presented. The theory was used to estimate 
heavy ion fragment momentum downshifts for relativistic oxygen and transverse 
momentum widths for relativistic lanthanum projectiles. The main new feature of 
this work was the identification of the imaginary component of the momentum 
transfer as the longitudinal collision momentum transfer. Finally, the use of the 
model as a mechanism for estimating collision impact parameters was described. 

The present theory is mainly applicable at intermediate or high energies 
because of the use of eikonal wavefunctions and the impulse approximation. At 
lower energies (below several hundred MeV/nucleon), the validity of straight line 
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trajectories and the assumption of a constant projectile velocity is questionable. 
Therefore, to compare theory with experiment at lower energies 22 revisions to the 
model are necessary. In particular, deceleration corrections to the constant velocity 
assumption are being developed. For incident energies greater than 1 AGeV, first- 
order deceleration corrections are small (< 1 percent). As the incident energy 
decreases, however, the first-order corrections increase significantly (over 50 percent 
at 100 A MeV), indicating that higher-order terms must be included. Work on this 
is in progress and will be reported when completed. 

The authors wish to thank Hank Crawford and Peter Lindstrom of Lawrence 
Berkeley Laboratory, and Frank Cudnotta of the Environmental Measurements 
Laboratory (USDOE) for useful comments and suggestions. This work was 
supported by NASA Grant Nos. NCCI-42 (F. K. and G. S. K.) and NAG-1 -797 


(J. W. N.) 
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Figure 1: Momentum transfer to the 16 0 projectile, as a function of impact 

parameter, for 2.1 AGeV oxygen colliding with a beryllium target. 


Figure 2: Target-averaged longitudinal momentum downshifts as a function of 

projectile fragment mass number for 2.1 AGeV 16 0 colliding with Be, C, 
Al, Cu, Ag, and Pb targets. The experimental data, taken from 
reference 2, are averaged over isotopes for each fragment mass. 


Figure 3: Transverse momentum widths as a function of fragment mass number 

for 1.2 AGeV 139 La colliding with a carbon target. The experimental data 
are taken from reference 16. 
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ABSTRACT 


Parameterizations of single nucleon emission from the electromagnetic interactions of cosmic 
rays with nuclei are presented. These parameterizations are based upon the most accurate 
theoretical calculations available today. When coupled with Strong interaction parameterizations, 
they should be very suitable for use in cosmic ray propagation through interstellar space, the 
Earth's atmosphere, lunar samples, meteorites and spacecraft walls. 



Introduction 


Galactic cosmic rays are very high energy particles confined to the region of our Milky Way 
galaxy. They consist of about 98% bare nuclei (stripped of all electrons) and about 2% electrons 
and protons (Simpson 1983). Of the nuclear component about 87% is hydrogen, about 12% is 
helium and the other 1% consists of heavier nuclei. Fe is the most abundant of these nuclei with a 
typical energy of about 1 GeV/N. Even though these heavy nuclei are not very abundant, they are 
very penetrating due to their large mass and high speed. 

An understanding of the interactions of galactic cosmic ray nuclei is important for several 
reasons: 

1) Knowledge of the cosmic ray spectrum at the top of the Earth's atmosphere and knowledge 
of the composition of the interstellar medium enables us to determine the cosmic ray spectrum at 
the source (Simpson 1983). 

2) Knowledge of the spectrum at the surface of the Earth and knowledge of the composition of 
the Earth's atmosphere enables us to determine the cosmic ray spectrum at the top of the 
atmosphere (Wilson, Townsend and Badavi 1987). 

3) The radiation environment inside a spacecraft, due to solar and galactic cosmic rays may be 
determined (Wilson and Townsend 1988). 

4) Studies of the history of extraterrestrial matter (such as lunar samples, meteorites and 
cosmic spherules and dust found in deep sea sediments) and also of the history of cosmic rays 
themselves can be made with knowledge of the production rate of various nuclides (Reedy 1987; 
Reedy, Arnold and Lai 1983). 

The basic nucleus-nucleus interaction that a cosmic ray undergoes can occur mainly via the 
Strong or Electromagnetic force. Strong interaction processes (Gyulassy 1981) have been studied 
extensively and quite recently the study of Electromagnetic processes in high energy nuclear 
collisions has begun (Bertulani and Baur 1988). 

In order to study the propagation of cosmic rays through interstellar space, the Earth's 
atmosphere or a spacecraft wall it is not enough to have only a good understanding of the nucleus- 



nucleus interaction mechanism. One must have an accurate theory of transport as well. Generally 
one uses a nucleus-nucleus interaction cross section as input to a transport computer code. These 
codes however can be very complex and therefore require simple expressions for the cross sections 
rather than the use of data bases or complicated theoretical expressions (Wilson and Townsend 
1988). Thus there has been a considerable effort to parameterize the cross section expressions so 
that the only required inputs are the nuclear energies and charge and mass numbers (Letaw, 
Silberberg and Tsao 1983; Silberberg and Tsao 1973; Townsend and Wilson 1986; Norbury, 
Cucinotta, Townsend and Wilson 1988; Wilson, Townsend and Badavi 1987). 

One approach to the parameterization of cross sections is to simply take all the available 
experimental data and fit a curve through it (Letaw, Silberberg and Tsao 1983; Silberberg and Tsao 
1973). Such an approach has certainly been useful and successful, but a much more satisfying 
parameterization would be one tied more directly to theory. It is the aim of the present work to 
obtain such a parameterization for the Electromagnetic (EM) part of the nucleus-nucleus interaction. 
One can then couple this with a similar theoretical parameterization of the Strong interaction 
process (Wilson, Townsend and Badavi 1987) to obtain a complete theoretical parameterization of 
the complete cross section. 

A preliminary parameterization of the EM process has already been presented (Norbury, 
Cucinotta, Townsend and Badavi 1988), which utilizes the Weizsacker-Williams (WW) method of 
virtual quanta (Bertulani and Baur 1988; Jackson 1975). However, since then the theory has been 
improved to include the effects of both electric dipole (El) and electric quadrupole (E2) interactions 
(Bertulani and Baur 1988; Norbury 1989a), which will henceforth be referred to as multipole 
theory in contrast to WW theory. These El and E2 effects modify the parameterization 
considerably. Also in the present work several different parameterizations are presented differing 
in degree of complexity. In addition much more data has become available with which to compare 
the parameterizations (Heckman and Lindstrom 1976; Olson, Berman, Greiner, Heckman, 
Lindstrom, Westfall and Crawford 1981; Mercier, Hill, Wohn, McCullough, Nieland, Winger, 
Howard, Renwick, Matheis and Smith 1986; Hill, Wohn, Winger and Smith 1988; Smith, Hill, 



Winger and Karol 1988; Hill, Wohn, Winger, Khayat, Leininger and Smith 1988; Hill, Wohn, 
Winger, Khayat, Mercier and Smith 1989; Norbury 1989b; Hill and Wohn 1989). The 
parameterizations to be presented below can then be combined with Strong interaction 
parameterizations such as the excellent parameterization by Wilson, Townsend and Badavi (1987). 
This combination should provide for much more accurate models of cosmic ray propagation 
through interstellar space, the Earth’s atmosphere and spacecraft walls. 

The present work will only consider single nucleon emission from cosmic ray nuclei. This has 
been shown to be the dominant electromagnetic process. Other particle emission processes such as 
two-neutron emission have much smaller probability (Hill, Wohn, Winger, Khayat, Mercier and 
Smith 1989), and will be studied in future work. 


ELECTROMAGNETIC THEORY 


The EM theory has already been discussed extensively (Bertulani and Baur 1988; Norbury 
1989a) and only a few relevant details will be given here. The total nucleus-nucleus EM cross 
section is written as 

a = gei + Oe2 

= J[Nei (E) oei (E) + N E2 (E) ct E 2 (E)] dE (1) 

where Nei (E) is the virtual photon spectrum (of energy E) of a particular multipolarity due to the 
projectile nucleus and Gei (E) + Ge 2 (E) is the photonuclear reaction cross section of the target 
nucleus. (In principle the above equation should include other EM multipoles, but their effect is 
much less important.) A less exact expression is given by WW theory as 

Gww (E) = jNww (E) [gei (E) + ge 2 (E)] dE (2) 

where Nww (E) is the WW virtual photon spectrum. Bertulani and Baur (1988) have shown 
that 

Nww (E) = N E i (E) = 1 J Z 2 a -i- Ko - 1 $ (Ki 2 - Kq 2 )] (3a) 

b Jl ^2 2 

N E 2(E)=l lz2 a-L [2 (1 - p 2 ) Ki 2 + £ (2 - (3 2 ) 2 KoK,-l (3 4 (Ki 2 - K 0 2 )] (3b) 

tit ^ i 


with 

_ Ebmin 

ypflic) 

where all of the Bessel functions K are functions of £. In the above equation E is the virtual 
photon energy, Z is the nuclear charge, a is the EM fine structure content, and b m in is the 
minimum impact parameter, below which the collision occurs via the Strong interaction. Also 



where c is the speed of light and v is the speed of the cosmic ray. 


(3 = ^ and 



The minimum impact parameter is given by 

bmin = Ro -1 (T) + Ro-i (P) - d (5a) 

where Ro-i are the 10 per cent charge density radii of the projectile and target and d is an adjustable 
overlap parameter. An excellent approximation to Ro-i is (Norbury, Cucinotta, Townsend and 
Badavi 1988). 

Ro-i = (1.18 AW + 0.75) fm (5b) 

where A is the nuclear mass number. 

Jackson has provided high and low virtual photon energy approximations as 

Nww(E) = li^aX[ln(l*i2i).lp 2 ] (6a) 

for small £, and 


Nww (E) - 1 Z 2 a (1 - 1 (3 2 ) exp (- 2$) 
E p2 2 


(6b) 


for large 

In equation (1) the El photonuclear cross section can be written in terms of the electric giant 
dipole resonance (GDR) cross section as 

o n 


Cei(E) = ■ 




1+[(E 2 -EJ dr ) 2 /E 2 Igdr] 

where Egdr is the energy of the peak in the GDR cross section, Tgdr is the width of GDR, and 


(7) 


0m = _2m_ (8) 

jJgdr/2 

with the Thomas-Reiche-Kuhn cross section (Levinger 1960) given by 

Otrk = MeV mb ^) 

where N and A are the neutron and mass numbers. The GDR energy is given by (Westfall, 
Wilson, Lindstrom, Crawford, Greiner and Heckman 1979) 


Egdr = tic [ — | (1 + u - 1 ^ e - 4 ' 3y - £ )1 

8J 1 + e + u 


- 1/2 


GO) 



with 


u=^A-W (11) 

and 

Ro = r 0 A 1/3 (12) 

where e = 0.0768, Q' = 17 MeV, J = 36.8 MeV, ro = 1.18 fm, and m* is 7/10 of the nucleon 
mass. Note that other expressions for Egdr such as 80 A' 1/3 (Bertulani and Baur 1988) provide 
very inaccurate results for light nuclei. Equation (10) is accurate for all mass regions. 

The E2 cross section is dominated by the giant quadrupole resonance (GQR). The main 
contribution to single nucleon emission (Bertulani and Baur 1988) comes from the isoscalar 
component given by (Bertrand 1976) 


<*E2 (E) = ■ 


Qewsr E 2 


1 + (E 2 - Eqqr) 2 / E 2 !' 2 gqr 

with the energy-weighted sum rule (EWSR) cross section 

0.22 7A m ^b MeV 4 
<*ewsr = f 

tcTgqr/2 

where f is the fractional exhaustion of the EWSR (Bertrand 1976) and 


(13a) 


(13b) 


EgQ r= aBt (i*) 

Finally, all of the above cross sections refer to total absorption cross sections. To obtain the 
reaction for proton or neutron emission they must be multiplied by the proton or neutron branching 
ratios. The proton branching ratio has been parameterized by Westfall et al (Westfall, Wilson, 
Lindstrom, Crawford, Greiner and Heckman 1979) as 

g p = Min [Z/A, 1.95 exp (-0.075 Z)] ( 15 a ) 

where Z is the number of protons and the minimum value of the two quantities in square brackets 
is to be taken. Assuming that only single nucleon emission occurs, the neutron branching ratio is 
gn = 1 - gp (15b) 



Comparison Between Theory and Experiment 


In the above paragraphs I have provided the basic equations to be used in the present work. 
However in analyzing the validity of the basic EM theory one uses only equations (3) - (5) and 
instead of equations (7) - (15) one uses actual experimental data for the photonuclear cross 
sections. A detailed study of the validity of this EM theory has been made (Norbury 1989a, b, c, 
d) and the results from this work are presented in Table 1, both for the WW theory and the separate 
El and E2 multipole theory calculations, and are compared to experimental data. A detailed 
discussion is to be found in Norbury 1989a, b, c, d, but the following features are to be noticed. 
Both WW and multipole theory give reasonably good results although multipole theory is 
somewhat better. It is found that electric quadrupole (E2) effects are not significant for proton and 
neutron emission from 12 C, ,6 0 or 18 0. However, E2 contributions are substantial for neutron 
emission from 59 Co, 89 Y and 197 Au, generally leading to improved agreement between theory and 
experiment Notable disagreements occur for 139 La projectiles (1.26 GeV/N) where the theoretical 
Oei + <*E2 are too big. Quadrupole effects improve the theoretical results for 16 0 projectiles at 60 
and 200 GeV/N, although the theoretical cross sections are still too small. In general it has been 
found (Norbury 1989a, d) that electric quadrupole effects are an important component in nucleus- 
nucleus collisions and that these effects can be calculated accurately. 


Parameterizations 


As mentioned above in testing the basic WW and multipole EM theory one uses experimental 
data for the photonuclear cross sections. However this is not a practical procedure for use in 
cosmic ray transport codes and instead my approach will be to use expressions (7) - (15). 

In the present work I shall discuss three separate parameterizations of the above EM theory for 
use in cosmic ray transport codes. These will be presented in decreasing order of accuracy, but the 
aim is to provide parameterizations that will be useful in different contexts. 

Parameterization #1 of Multipole Theory 
This is the most accurate parameterization and uses the following equations: 

1) Equation (1) is used for the total nucleus-nucleus EM cross section. The integration is 
done numerically using the Trapezoidal Rule. 

2) Equations (3) are used for the virtual photon spectra Nei (E) and Ne 2 (E). 

3) Equations (5) are used for the minimum impact parameter with the overlap parameter 
adjusted to give the best fit to data at d = -1.5 fm. 

4) Equations (7) - (14) are used for the photonuclear cross sections. 

5) The width Tgdr in equations (7) and (8) is set at 

r G DR= lOMeV for A <50 

= 4.5 MeV for A > 50 (16) 

and Tgqr in equations (13) is set at 

Tgqr = 2.5 MeV for A > 180 

= 4.5 MeV for 70<A<180 
= 5.5 MeV for 19 < A <70 

= 3.0 MeV for A < 19 (17) 

These values for Tgrd are discussed in Norbury, Cucinotta, Townsend and Badavi 
(1988) and for Tgqr in Bertrand (1976). 

6) The fractional exhaustion of the Energy-Weighted Sum Rule in equation (13b) is given 
by (Bertrand 1976) 

f = 0.9 for A> 100 
= 0.6 for 40 < A < 100 
= 0.3 for 40 < A 

7) The proton and neutron branching ratios are given by equations (15). 


( 18 ) 



The results of the above parameterizations are given in Table 1. It can be seen that it agrees 
extremely well with the multipole theory. Thus I regard this parameterization #1 of the multipole 
theory as describing very accurately the most advanced state-of-the-art EM theory. Agreement 
between this parameterization and experiment is, of course, of the same quality as between the 
multipole theory and experiment. 

Parameterization #1 of WW Theory 

WW theory gives a simpler treatment of the virtual photon field and is included here for the 
sake of completeness. The only difference between parameterization #1 of WW theory and 
parameterization #1 of multipole theory is that equation (2) is used for the total cross section 
instead of equation (1). Results are listed in Table 1 and are fairly comparable to the 
parameterization #1 of the multipole theory. 

Parameterization #2 of Multipole Theory 

A difficulty that might occur in some cosmic ray transport theories is the necessity of having to 
do a numerical integration in equation (1) every time a is to be evaluated. To get around this, 
parameterization #2 is based on the technique of Bertulani and Baur (1988). This involves taking 
NEi (E) outside of the integral in equation (1) and evaluating Nri (E) at Egdr (see equation 10) and 
Ne 2 (E) at Eqqr (equation 14). The remaining integral is evaluated from sum rules. That is 
(Bertulani and Baur 1988), equation (1) becomes 

CJ - Nei (Egdr) 1 <J6i (E) dE + N E2 (Egqr) e£q R J^^dE (19) 

E z 

with the sum rules 

J cei (E) dE - 6Q-HZ. m b MeV < 20a ) 

A 

and 

J Pe 2 (E) ^ =f Q.22 ZA 2 ^ mb MeV' 1 (20b) 

1000 

Bertulani and Baur (1988) claim that this is an accurate procedure. However, I found it necessary 
to change d to d = -2.4 fm (see equation 5a) in order to give good comparison to experiment. 


In the present parameterization #2 of multipole theory items 1) - 3) of parameterization #1 were 
changed to those discussed in the preceding paragraph. Note especially that a numerical integration 
is no longer necessary. Items 4) - 5) are no longer relevant. Items 6) - 7) remained the same. 
Results are again listed in Table 1. With the new value of d = -2.4 fm parameterization #2 agrees 
well with parameterization #1 (which used d = -1.5 fm). 

Parameterization #2 of WW Theory 

WW theory is again included for completeness. In this case equation (2) was replaced with 

Cww * Nww (Egdr) J Cei (E) dE + Nww (Egqr) Egq R (21) 

with the same sum rules in equations (20). Results are listed in Table 1. 

Parameterization #3 

Parameterizations #1 and #2 require the evaluation of Bessel functions as indicated in equations 
(3) for the virtual photon spectra. In the interest of providing an even simpler parameterization that 
could be used on a pocket calculator for rough estimates of the cross section, a third 
parameterization is presented. The E2 cross section was ignored and equations (1) or (2) were 
replaced with 

a = Nww (Egdr) J Ofei (E) dE (22) 

Note that this is identical to neglecting the GQR in equation (1). The sum rule in equation (20a) 
was used for the integral. Nww (Egdr) was evaluated using equations (6), with (6a) used for 
t, < 0.5 and (6b) for £, > 0.5. This prescription avoids the evaluation of Bessel functions and 
almost allows one to calculate a in one's head. In this case the value of d was d = +1.0 fm. Items 
4) - 6) are not relevant and item 7) was again used. The results are presented in Table 1 and are 
seen to give surprisingly similar results to the other parameterizations. 




As discussed in previous work (Norbury 1989a, d) the multipole theory is generally more 
accurate than WW theory. This is also true for the above parametrizations as can be seen from 
Table 1. 

However, WW theory and multipole theory do not describe 18 0 very well, and the 
parameterizations are even worse. I trace this to the fact that the branching ratio equations (15) do 
not work well for nuclei off the stability curve. 

Both WW theory and multipole theory do not describe 197 Au very well either, but the 
parameterizations do a somewhat better job due to the choice of the overlap parameter d. There 
seems to be a problem also for very high energies especially 200 GeV/N. 

Apart from these problems the multipole theory and multipole parameterizations (#1, #2 and 
#3) seem to describe the data quite accurately. 

As regards which parameterization to use, they all seem to do an equivalent job in describing 
the data. This of course is because a different value for d was chosen for each. Even the 
parameterization #3 does quite well, although it is a little high for nucleon emission from the lighter 
nuclei. 

Given the above problems with 18 0, 197 Au and 200 GeV/N I recommend that the above 
parameterizations be used 0 only with nuclei on the stability curve, *0 for nuclei lighter than 197 Au 
and '“) for energies less than 10 GeV/N. These requirements should not be too restrictive in 
Cosmic Ray work because most nuclei have energies of around 1 GeV/N and the most abundant 
nuclei are not much heavier that 56 Fe (Simpson 1983). Having to deal only with nuclei on the 
stability curve is probably the most severe restriction. 

Parameterizations #1, #2, #3 decrease in order of accuracy, but, as discussed above, not by 
very much. I would recommend using the most accurate parameterization (#1), but if one's 
computer codes are such that it would save CPU time by using either #2 or #3, then I would 
recommend their use. However, one should perhaps be careful about using parameterization #3 


for light nuclei. I recommend the multipole parameterizations, but I do not recommend the use of 
the WW parameterizations. 

Finally, by combining the above EM parameterizations with the Strong Interaction 
parameterization of Wilson, Townsend and Badavi (1987), which is not subject to the same 
restrictions as above, transport of cosmic rays through matter can be described very accurately. 
Future work will involve parameterization of both multiple nucleon emission (a much smaller 
effect) and also neutron branching ratios for nuclei off the stability curve. 
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Table 1 Theoretical (Norbury 1989 a, c, d), Parameterized (this work) and Experimental (Heckman and Lindstrom 1976, Olson et al. 1981, Mercier et al. 1986, 
Smith et al. 1988, Hill et al. 1988, 1989) EM Cross Sections, all in units of milli-barn (mb). 
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Abstract 


Calculations are presented for electric quadrupole excitations in relativistic nucleus- 
nucleus collisions. The theoretical results are compared to an extensive data set and it is 
found that electric quadrupole effects provide substantial corrections to cross sections, 
especially for heavier nuclei. 


L INTRODUCTION 


The search for a fundamentally new state of matter in the form of a Quark-Gluon 
Plasma *) has stimulated the production of very high energy nuclear beams. The hope is to 
observe the Quark-Gluon Plasma in a relativistic nucleus-nucleus collision. At the 
Berkeley Bevalac a variety of light nuclei such as 12 C, 1 ^0 and 20 Ne can be accelerated up 
to energies of 2. 1 GeV/N and heavier nuclei such as 139 La and 238 U can be accelerated to 
1.26 and 0.96 GeV/N respectively. At Brookhaven, New York, 16 0 beams are available at 
14.6 GeV/N and at the CERN SPS, in Geneva, beams of 1 ^o and 32 S are both produced at 
60 and 200 GeV/N. The Relativistic Heavy Ion Collider (RHIC) is expected to produce 
two colliding beams at 100 GeV/N to give a total center-of-mass energy of 200 GeV/N, 
which corresponds to a single beam energy of 21 TeV/N. Grabiak 2 ) has pointed out that 
nuclear beams of 3.5 TeV/N and 8 TeV/N may be possible at the CERN Large Hadron 
Collider (LHC) or the Superconducting Super Collider (SSC). By way of comparison, the 
majority of Galactic Cosmic Rays have energies 3 ) of about 1 GeV/N, with a range 
typically from 10 MeV/N to 1 TeV/N. However, the JACEE collaboration 4 ) has made 
observations as high as 1000 TeV/N. 

Nucleus-nucleus reactions proceed mainly through either the Strong or Electromagnetic 
(EM) interactions. Historically, Strong interaction processes have been the main object of 
study 5 ), however with the availability of the above high energy nuclear beams there has 
been a resurgence of interest in EM interactions in relativistic nucleus-nucleus collisions. 6) 

The primary theoretical tool for studying these relativistic EM processes has been via 
the Weizsacker-Williams (WW) method 6 ‘ 7 ) of virtual quanta. The nucleus-nucleus total 
EM reaction cross section is 

Ct = j Nww (Ey) tf(Ey) dEy (1) 

where Ey is the virtual photon energy, N^y (Ey) is the WW virtual photon spectrum and 
c(Ey) is the photonuclear reaction cross section. For high accuracy it is important to use 
experimental photonuclear data for o(Ey). (For an excellent compilation of photoneutron 


data see reference 8.) However, a more exact formulation of a involves a breakdown into 
the various EM multipolarities such as electric dipole (El), electric quadrupole (E2), 
magnetic dipole (Ml) etc. The most important contributions to a are from El and E2 so 

that 


a - a m + a E2 


I 


[NEi(Ey)OEi(Ey) + N E2 (E y )a E 2(E y )] dEy 


( 2 ) 


where Ngi (Ey) is the virtual photon spectrum of a particular multipolarity due to the 
projectile nucleus and Ogi (Ey) is the photoriuclear reaction cross section of the target 
nucleus. Bertulani and Baur 6) have derived expressions for Ngi (Ey) and found that the 
electric dipole spectrum is the same as the WW spectrum, i.e. Ngi (Ey) = N w (Ey). 
Furthermore at very high projectile energies all Ngi (Ey) and Nmi (Ey) are equal so that 
equation (1) is seen to be a very high energy approximation to all multipolarities included in 
equation (2). Bertulani and Baur 6 ) have made a crude estimate of the EM cross section 
using equation (2) but they pulled Nei (Ey) and Ne 2 (Ey) outside the integral and evaluated 

them at a single energy and used sum rules to evaluate J °Ei (Ey) dEy. However, this 

procedure did not yield very accurate results. Thus I undertook a more exact study 9 ) 
leaving equation (2) as it stands, and using experimental data for the photonuclear cross 
sections by defining 

(3) 

^E1 (Ey) = <3 eX pt. (Ey) — <TE2(Ey) 

where a ex pt. (Ey) is the experimentally measured photonuclear cross section and (Te 2 (Ey) is 
a theoretical calculation based on a Lorentzian shape for the Electric Giant Quadrupole 
Resonance (GQR). Details for this procedure can be found in reference 9. As was noted 
in that reference, the above procedure yields very accurate values for the sum agi + cj£2 
(which is to be compared to nucleus-nucleus reaction experiments) even though the GQR 
parameters are uncertain. The basic reason for this, as can be seen from equation (3), is 


that an under (over) estimate in Oe 2 (Ey) will give an over (under) estimate in agi (Ey), so 
that the combined <7 ei + <5E2 in equation (2) will not change very much. 

In reference 9 a detailed study of El and E2 was undertaken for the reaction 89 Y (RHI, 
X) 88 Y where RHI refers to various Relativistic Heavy Ions and X is anything. It was 
found that E2 effects account for a considerable fraction of the cross section, and that 
inclusion of E2 (via equation 2) provides improved agreement with experiment over the 
WW method. Given this situation, it was decided to compare this theoretical approach to 
as much experimental data as possible. Thus the present work involves a comparison to 
neutron emission from 89 Y, 197 Au, 59 Co and neutron and proton emission from 12 C, 16 0 
and 18 0 which includes both electric dipole and quadrupole effects. This complements 
earlier work 7 ) which involved an extensive comparison of WW theory to experiment. 
CALCULATIQNAL METHOD 

The basic calculational method is outlined in reference 9 and the discussion will not be 
repeated here. Also reference 7 includes a very detailed summary of which photonuclear 
data were used for a expt . (Ey) in equation (3). The same data is used in the present work. 
All isoscalar GQR parameters were taken from the compilation of Bertrand 10 ) and are 
listed in Table 1. As mentioned in the Introduction, even though these parameters are 
somewhat uncertain the total EM cross section ogi + 0£2 is expected to be very accurate 9 ) 
due to the subtraction procedure of equation (3). The most inaccurate results would be 
expected for the 12 C, 16 0, 18 0 GQR parameters where the isoscalar GQR is fragmented 
into several components 10). Only a single Lorentzian 9 ) was 



Table 1: Isoscalar Giant Quadrupole Resonance (GQR) Parameters taken from the 
compilation of Bertrand 10 \ E is the GQR resonance excitation energy, T is the full-width 
at half maximum and f is the fractional depletion of the Energy Weighted Sum Rule. (The 
GQR of light nuclei are fragmented into several peaks, so that the parameters below 
represent an estimated average value.) 


Nucleus 

E 

(MeV) 

r 

(MeV) 

f 

12C 

27. 5 a 

3.0 b 

0.4 d 

16o 

22.0 b 

3.0 C 

0.4 d 

18o 

24.0 3 

3.0 C 

0.4c 

59Co 

16.3 b 

5.6 b 

0.61 b 

89y 

13. 8 b 

3.2 b 

0.55 d 

192 Au 

10.8 a 

2.9 b 

0.95 d 


a E is calculated from 63 A' 1 / 3 
b best value taken from Table 4 of ref. 10 
c estimate 

d from Fig 23 of ref. 10 

used in the present work. However oe 2 is found to be quite small for these nuclei (see 
below) so that my conclusion that the calculated cjei + OE2 is accurate remains valid. 

For the nuclei 12 C, 16 0 and 18 0, proton (p) emission occurs as well as neutron (n) 
emission. Thus equation (3) needs to be modified to incorporate the branching ratio. I 
assume that the excited nucleus decays only by proton or neutron emission and that the 
(photon) energy dependent neutron branching ratio is defined as 

^expL(Ey>n) + ^expLCEyiP) 


so that 


(5) 


CE2 (Ey, n) = f n (Ey) OE2 (Ey) 

where Ge2 (Ey) is the photonuclear GQR cross section. Thus for proton and neutron 
emission equation (3) becomes 

Cei (Ey, n) = Gexpt. (Ey, n) - f n (Ey) ge2 (Ey) (6a) 

and 

Oei (Ey, p) = G exp t. (Ey, p) - [1 - f n (Ey)] Ge2 (Ey) (6b) 

Equations (4) - (6) were used for nucleon emission from 12 C, l6 0 and 18 0. For 59 Co, the 
(y, p) cross section is not available and so a constant value of f n = 0.7 (suggested from 
reference 11) was used. For 89 Y and 197 Au I used f n = 1.0. 

RESULTS AND DISCUSSION 

The calculated results are listed in Table 2, along with the experimental results of 
various groups. 12-16) <j E1 + <j E2 is the calculated result to be compared with the data 
Gexpt- Also listed are the results of WW calculations. 7 ) In all cases two theoretical cross 
sections are listed. The first is calculated using an expression for the minimum impact 
parameter as 

b m in = Ro i(T) + R 0 i (P) (7) 

where Roi represents the 10-percent charge density radius 7 ) of the target or projectile. 

The second theoretical cross section listed in parentheses in Table 2 uses b m j n given by Hill 
et al. 14-16) ^ 

bmin = r 0 [A p V3 + A T W - X(A p l/3 + Ar 1 ' 3 ] (8) 

where r 0 = 1.34 fm and X = 0.75. (Note that my WW calculations disagree with earlier 
results of Hill et al. 14-16) due to an error in their calculations. 19-20)) 

There are several features readily apparent from Table 2; 

(0 Gei + Ge 2 is always larger than Gww- However, for nucleon emission from 12 C, 
16 0 and 18 0 this difference is never larger than about 4%, but for neutron emission from 
59 Co, 89 Y and 197 Au the difference is much larger varying between about 7% - 15%. 



“) For nucleon emission from 12 C and 16 0 both oei + <*E2 and agree with 
experiment for both choices of bm^. 

“0 For nucleon emission from 18 0 both ctei + GE2 and disagree with experiment 
for both choices of bmin. y actually gives slightly better agreement but not by a 
significant amount. 

iv ) For neutron emission from 197 Au, Gei + <*E2 is significantly closer to experimental 
values than is Oww> although for most cases it still lies outside the error bars. An 
exception however is a much poorer agreement for 139 La (see also references 19 and 20). 
Significant discrepancies with 197 Au data have been noted previously for WW theory. 7 ) 

v ) For neutron emission from 89 Y, oei + &E2 is in much better agreement with 
experiment than is oww- This is especially true for the 40 Ar and 56 Fe projectiles. 

vi ) For 59 Co, oei + <TE2 is again better for 20 Ne, although slighdy worse for ^Fe. As 
above the agreement for the 139 La projectile is significantly poorer. 

S U MMARY AND CO NC LUSIO NS 

Calculations have been made for nucleon emission via EM dissociation in relativistic 
nucleus-nucleus collisions. Results are presented for Weizsacker-Williams theory and also 
for separate electric dipole and quadrupole components. The theories have been compared 
to an extensive data set. It is found that electric quadrupole (E2) effects are not significant 
for proton and neutron emission from 12 C, 16 0 or 18 0. However, E2 contributions are 
substantial for neutron emission from 39 Co, 89 Y and 197 Au, generally leading to improved 
agreement between theory and experiment. Notable disagreements occur for 139 La 
projectiles (1.26 GeV/N) where the theoretical oei + &E2 are too big. Quadrupole effects 
improve the theoretical results for 16 0 projectiles at 60 and 200 GeV/N, although the 
theoretical cross sections are still too small. 

In general it has been found that electric quadrupole effects are an important component 
in nucleus-nucleus collisions and that these effects can be calculated accurately. 
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Table 2 Calculated results, cei + Oe 2 and cww> compared to experiment. 12-16) t wo theoretical cross sections are listed, 
first set uses b m i n given by equation (7) and the second set (in parentheses) uses b m i„ given by equation (8). 



15n I 30±16 | 42 (46) | 42 (46) | 1 (2) | 43 (48) 
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1.05 U C -2±5 0.3 (0.4) 0.3 (0.4) 0(0) 0.3 (0.4) 

" n B -2±5 0.5 (0.6) 0.5 (0.6) 0(0) 0.5 (0.6) 
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Abstract 


The first complete calculations of electric quadrupole excitations in relativistic nucleus- 
nucleus collisions are presented herein. Neutron emission from 89 Y is studied and 
quadrupole effects are found to be a significant fraction of the cross section. 


PACS 25.70.Np 



Nucleus-nucleus collisions proceed predominantly via the Strong and 
Electromagnetic (EM) forces, both of which have been studied extensively. 1 * 2 ) The EM 
interaction consists of many multipoles such as electric dipole (El), electric quadrupole 
(E2), magnetic dipole (Ml) etc. The electric dipole is the most important of these and this 
has been the only EM multipole for which calculations have been made and compared to 
experiment in relativistic nucleus-nucleus collisions. In this paper I present the very first 
accurate calculations of the electric quadrupole effect. 

In nucleus-nucleus collisions the Strong interaction dominates the cross section at 
impact parameters approximately less than the sum of the nuclear radii, ie. for impact 
parameters smaller than 

bmin ~ Ro-l(T) + Ro-i(P) (1) 

where Ro i represents the 10-percent charge density radius 3 * 4 ) of the target or projectile. 
(Other expressions for b m in are possible, 2 > 5 ) but for the sake of simplicity they are not 
discussed here. For impact parameters larger than b m in> the interaction occurs via the EM 
force and the EM cross section is calculated via 


° = °El + a E2 


= J[N e 1 (E) a m (E) + N e 2 (E) g £2 (E)] dE 


( 2 ) 


El v ' El v ' E2' ' E2 V 

where Ngi (E) is the virtual photon spectrum (of energy E) of a particular multipolarity due 
to the projectile nucleus and OEi (E) is the photonudear reaction cross section of the target 
nucleus. (In principle the above equation should include other EM multipoles, but their 
effect is much less important.) 

All previous comparisons between theory and experiment 5 ' 9 ) have only included 
the electric dipole effect using Nei(E) calculated from Weizsacker- 


Williams (WW) theory 2 > 10 ) where 

j N WW (3) 

where o(E) is the experimentally measured photonuclear reaction cross section. Nww (E) 
is equal to Nei (E) 2 ) so that WW theory does not include the quadrupole component 

There are difficulties in evaluating a in equation (2). To be as accurate as possible, 
one should use experimental values for OEi (E) and OE2 (E), fold them into the energy 
dependent spectra Nei (E) and Ne 2 (E) and integrate the whole expression numerically. 
Bertulani and Baur 2 ) have made a crude estimate of the EM cross section using equation 
(2). However they pulled Nei and Ne 2 outside of the integral and evaluated them only at a 
single energy corresponding to a theortical estimate of the peak in the El and E2 cross 
section. The remaining J o(E) dE for El and E2 were evaluated using theoretical sum 

rules. This procedure led, for example, to a total cross section a of 839 mb and 266 mb 
for the reactions 197 Au ( 56 Fe, X) 196 Au and 89 Y ( 56 Fe, X) 88 Y respectively (at 1.7 
GeV/N), whereas the measured cross sections are 601 ± 54 mb and 217 ± 20 mb 
respectively 5 X Given such a discrepancy, I decided to retain the energy dependence in 
N(E) by doing a numerical integration, as described above, using, where possible 
experimental photonuclear cross sections, as detailed below. 

In the present work, results are presented for the reaction 89 Y (Projectile, X) 88 Y 
only; the major point being simply to illustrate the importance of E2 effects using an 
accurate calculation. Results for other nuclei such as 12 C, 16 0, 18 0, 59 Co and 197 Au and 
detailed comparisons to data will be presented elsewhere. 

For best accuracy I have followed the suggestions of Berman et al. concerning 
which photonuetron reaction data to use for 89 Y. Following their suggestion I have used 
the Saclay 12 ) data but multiplied by a factor of 0.82. (The data actually stops at 27 MeV 
and a smooth extrapolation was used to estimate the small amount of remaining data 
beyond this energy.) However, all experimental photoneutron data consists of oei (E) 


plus oe 2 (E) and a way must be found to separate out these components so that they can be 
inserted into equation (2) and numerically integrated. This separation was achieved by 
using a theoretical calculation 13 ) of the isoscalar component of the electric giant 
quadrupole resonance (GQR) 




a E 
EWSR 


'^-EGQR^r 2 


(4) 


with the energy- weighted sum rule cross section 


a EWSR f 


0.22 Z A^jib MeV 1 
ttr/2 


(5) 


The parameters in the above expressions were taken from Bertrand. 14 ) For 89 Y the width 
T is 3.2 MeV, the energy of the giant quadrupole resonance Egqr is 13.8 MeV and the 
fractional exhaustion of the EWSR f is 55%. (Note that the 89 Y nucleus is approximately 
spherical, thus justifying the use of a single Lorenman in equation (4).) The expressions 
for oe 2 (E) in equations (4) and (5) were used in equation (2). The dipole cross section 
was determined by subtracting CE2 (E), as given above, from the experimental cross 
section <J ex pt. (E) of Lepretre 12 ) as in 

CfEl (E) = OexpL (E) - a E2 (E) (6) 

where a eX pt. is 0.82 times the Lepretre 12 ) cross section. Then oe 2 (E) and oei (E) were 
inserted into equation (2) and Nei (E) and Ne 2 (E) were taken from expressions derived by 
Bertulani and Baur 2 ) The integrals in equation (2) were performed numerically to give the 
EM nucleus-nucleus cross sections. Because of the use of equation (6) , uncertainties in 
the GQR parameters (even if they were as large as ± 2MeV in T and Egqr and ± 20% in f) 
do not change the total calculated EM cross section oei + OE2 (which is compared to data 
in Table 1) by more than 4%. Thus the calculations presented herein are expected to be 
very accurate even if the quadrupole parameters are uncertain. 



Results for the reaction 89 Y (projectile, X) 88 Y are presented in Table 1 and 
compared to the experimental measurements of Mercier et al. 5 ) Both individual dipole ogi 
and quadrupole ag 2 cross sections are presented as well as their sum a which is to be 
compared to the data. Also presented are results obtained using WW theory. (Note that the 
EM calculations using WW theory in reference 5 are not correct. 15, 16)) in all calculations 
bmin from equation (1) are used. The 10-percent charge radii 3 > 4 ) are also listed in Table 1. 

One can see that WW theory agrees with experiment for the 12 c and 20Ne 
projectiles and is reasonably close for the 40 Ar and 56pe projectiles. Agreement could be 
reached by using a different expression for b m in. However a detailed study of treating 
bmin as an adjustable parameter will be reported elsewhere I 7 ). The quadrupole cross 
sections 0 E 2 are all seen to be about 10% of the dipole cross sections ctei, and for all 
reactions ogi + ctE2» is about 7% bigger than <7ww- This is because Ng2 (E) is always 
larger than Nei(E) so that the quadrupole photonuclear component oe 2 (E) is enhanced 
over the dipole. Adding the quadrupole now gives improved agreement between theory 
and experiment. 

In summary, the first accurate calculations of electric quadrupole effects in 
relativistic nucleus-nucleus collisions are reported. For the reaction 89 Y (Projectile, X) 88 Y 
the quadrupole cross section is about 10% of the dipole cross section, and thus I conclude 
that electric quadrupole effects are an important consideration in the analysis of nucleus- 
nucleus collisions. 
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Table 1 EM Cross Sections for the Reaction 89 Y (Projectile, X) 


The Photoneutron Cross Section measured by Lepretre et al 12 ) (but multiplied ^ by 0.82) 
was used in the analysis described in the text. The 10 percent charge radius used for 89 Y is 
6.02 fm and the GQR parameters (see text) arc f = 0.55, T = 3.2 MeV, Eqqr = 13.8 MeV. 
The 89 Y (y, n) threshold is at 1 1.0 MeV. Calculations are made for Weizacker- Williams 
theory (oww) and also individual El & E2 multipole cross sections are calculated. The 
total cross section ctei + <J£2 is to be compared to experiment. All calculations use the 
minimum impact parameter given by b m i n = R o-l (P) + R 01 (X). 


Projectile 

Ro-l (P) 

(fm) 

Energy 

(GeV/N) 

°expt. ( m b) 
(ref. 5) 

(mb) 

<?E1 

(mb) 

<*E2 

(mb) 

ctei + OE2 
(mb) 

12 C 

3.30 

2.1 

9± 12 

12 

12 

1 

13 

20 Ne 

4.00 

2.1 

43 ±12 

32 

31 

3 

34 

40 Ar 

4.72 

1.8 

132 ±17 

90 

88 

9 

97 

56 Fe 

5.24 

1.7 

217 ± 20 

175 

171 

16 

187 
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Abstract 

Previous analyses of the comparison of Weizsacker-Williams (WW) theory to experiment 
for nucleon emission via electromagnetic (EM) excitations in nucleus-nucleus collisions 
have not been definitive because of different assumptions concerning the value of the 
minimum impact parameter. This situation is corrected by providing criteria that allow one 
to make definitive statements concerning agreement or disagreement between WW theory 
and experiment. 

PACS: 25.70. N p 



1. INTRODUCTION 


Collisions between relativistic nuclei can occur via the Strong or Electromagnetic 
interaction. There is an enormous literature on processes induced by the Strong force 1 ), 
but relatively few studies have been carried out on the Electromagnetic (EM) aspects of 
relativistic nucleus-nucleus collisions. 2 ' 25 ) This situation is surprising given the richness 
of applications of EM effects. These effects are of importance for the following reasons: i) 
EM interactions between relativistic nuclei are interesting in their own right; ii) they will 
form a significant background to the formation of a quark-gluon plasma at ultrarelativistic 
energies; iii) other applications in physics such as subthreshold pion production 3 ) iv) 
astrophysical applications 4 ); v) interference effects between Strong and EM amplitudes 5 ); 
vi) studies of virtual photon theory 4 > 6-10 ) and vii) applications in space radiation effects **) 
Bertulani and Baur 25 ) have written an outstanding review article on EM effects in nucleus- 
nucleus collisions to which the reader is referred. 

The present paper is concerned with nucleon emission via electromagnetic 
dissociation in relativistic nucleus-nucleus collisions. The first experiment of this kind 
was performed by Heckman and Lindstrom 15 ) looking at excitations in 12 C and 16 0 
projectiles at energies of 1.05 and 2.1 GeV/N on a variety of targets ( 12 C, 27 A1, ^Cu, 

108 Ag, 208 Pb). Measured EM cross sections for nucleon emission ranged from 0 to 
50 mb. Olson et al 13 ) later measured excitation of 18 0 projectiles at 1.7 GeV/N on 48 Ti, 
208 Pb and 238 U with cross sections up to 140 mb. Studies of ^ 7 Au and 59 Co target 
excitation 14 > 20-22) were j ater reported with cross sections all the way up to 1970 mb for 
139 La projectiles at 1.26 GeV/N. Lighter projectiles were also used 14 > 20 * 23 ) with smaller 
cross sections. All studies mentioned so far have been for projectile energies less than or 
equal to 2. 1 GeV/N. The measurements were made at the Berkeley Bevalac. However, 
some very interesting measurements have also been made for 197 Au target excitation at the 
CERN SPS using projectiles at 60 and 200 GeV/N with cross sections of 820 and 
440 mb respectively. All the above data is summarized in Tables 1 and 3. 
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The authors of the above experiments have generally made a comparison of their 
data to theoretical predictions based on the Weizacker-Williams (WW) method of virtual 
quanta. 6 ’ 12 » 15 > 25 ) The basic idea is that the vertical photon spectrum N(E) of one nucleus 
is calculated from WW theory. 6 ) This is folded into the photonuclear cross section a v (E) 
for processes induced in the excited nucleus and then numerically integrated over energy to 
give the total EM nucleus-nucleus cross section 

<^ M = |o v (E)N(E)dE (1) 

Expressions for N(E) are given by Jackson 6 ) for the WW theory. These expressions 
include a minimum impact parameter b, below which EM interactions are supposed not to 
take place; the interaction proceeding via the much stronger nuclear force. One might 
naively expect b to be just the sum of the two nuclear radii. Note also that in WW theory 
N(E) is the same for all EM multipoles 25 ) and so Ov(E) does not need to be divided into its 
constituent multipoles. It is partly for this reason that the present paper concerns itself 
exclusively with WW theory. Alternative theories for N(E) 7 > 8 - 25 ) require that Oy(E) be 
divided into its constituent multipoles which is very involved and beyond the scope of the 
present paper, although work is proceeding in this direction. Here I wish to analyze WW 
theory only. 

Agreement between WW theory and experiment has generally been claimed to be 
good 1 15, 20-23, 25) an( i U pon reading the literature on the subject one is left with the 
impression that WW theory is an accurate theoretical tool. However, if one examines the 
calculations and data more carefully one is lead into some serious doubts. The following 
concerns arise. 

!) Table 1 provides a list of some experimental cross sections versus the theoretical 
calculations presented by various authors. As can be seen there is, in fact, very noticeable 
disagreement between experiment and WW theory, which might lead one to conclude that 
WW theory is no good at all. 
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2 ) Each research group 13 " 15 - 20 * 23 ) uses a different procedure for determining the 
minimum impact parameter, and hence agreement between theory and experiment depends 
to some extent on whose impact parameter one chooses. The differences in this parameter 
are in fact large enough to change agreement between theory and experiment. 

3 ) The calculation of the EM nucleus-nucleus cross section depends heavily on the 
photonuclear cross section Cv(E) used in equation (1). These authors have always used 
experimental data for the photonuclear cross section o v (E). The trouble is that various 
experimental data for Cv(E) for a particular nucleus are often in disagreement 26 ) and this 
can lead to significant differences in the calculated nucleus-nucleus EM cross section 
depending on whose data one chooses for Ov(E). 

For the above reasons it was decided to re-analyze the WW method applying a 
single method to all existing data paying particular attention to the following points. 

A ) Some recent articles have made a detailed study of the conflicting experimental 
photonuclear data 27 > 28 ) and have made recommendations concerning what is the correct 
data. The a v (E) used in the present work is that recommended by these studies 27 - 28 ) for 
12 C, 16 0 and 197 Au. For 18 0 and 59 Co the data of references 33 and 32 are expected to be 
very accurate. The two sets of conflicting data 30, 3 1) f or 89y are both used herein for two 
separate nucleus-nucleus calculations. 

B ) Some incorrect calculations have been presented in the literature. 2 ®' 22 ) They are 
listed correctly herein. (See references 37, 38 for a discussion of these corrections). 

C) To get around the problem of different possible choices of the minimum impact 
parameter b, it was decided to use it simply as an adjustable parameter fitted so that the 
theoretical WW cross section is equal to the experimental one. If b is a reasonable 
(unreasonable) value then one can definitely say that WW theory does (does not) agree with 
experiment. A reasonable value would be the sum of the two nuclear radii, whereas an 
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Table 1 Selected Data and Theoretical Ca lculations Showing Marked 
Disagreement between Theory and Experiment 


Projectile 

Target 

Energy 

(GeV/N) 

Final 

State 

^expt. 

(mb) 

from Reference 

Reference 

^theory 

(mb) 

from Reference 

18 0 

Ti 

1.7 

170 

8.7 ± 2.7 

13 

12.5 

" 

it 

tt 

17 N 

-0.5 ± 1.0 

" 

2.4 

n 

Pb 

n 

16o 

65.2 ± 2.3 

it 

55.2 

ii 

U 

ii 

170 

140.8 ±4.1 

tt 

167 

ii 

M 

n 

17 N 

25 1 ± 1.6 

it 

29.2 

n 

II 

it 

160 

74 .3 ± 1.7 

ii 

68.1 

12C 

192 Au 

2.1 

196Au 

75 ± 14 

20 

45 

20Ne 

n 

ii 

ii 

153 ± 18 

n 

121 

139La 

it 

1.26 

M 

1970 ± 130 

22 

2340 

56 Fe 

89y 

1.7 

88y 

217 ± 20 

20 

248 

56 Fe 

59Co 

ti 

58 Co 

88 ± 14 

it 

122 

139La 

tt 

n 

tt 

280 ±40 

22 

430 

16o 

192 Au 

60 

196Au 

280 ± 30 

21 

220 

n 

tt 

200 

ii 

440 ± 40 

n 

300 

20Ne 

tt 

2.1 

195Au 

49 ± 15 

23 

14 

40Ar 

tt 

1.8 

II 

76 ± 18 

tt 

38 

139La 

II 

1.26 

It 

73 ± 13 

it 

238 
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unreasonable value might be close to 0 or very much larger than the sum of the radii. Note 
that there will be cases where one cannot make a definite statement concerning agreement 
between theory and experiment. However, the advantage of the above criterion is that it 
does provide for an unambiguous comparison between theory and experiment. One will 
not be left wondering whether the agreement or disagreement with experiment is due to a 
choice of parameters. The value of the fitted parameter b will provide either a definitive 
conclusion regarding agreement or disagreement and in the cases in which such a 
conclusion is not possible it will be clear why 'such a conclusion is, in fact, not possible. 

In summary, the present work goes beyond other studies 13-15, 20-23, 25) j n 4 
respects. First, it provides a comprehensive comparison to all existing nucleon emission 
data. Second, problems due to inaccurate photonuclear input data are avoided. Third, 
previous incorrect calculations are corrected. Fourth, obtaining the value of b needed to fit 
the data enables one to make definitive statements concerning the agreement with WW 
theory, thus sidestepping the problem of comparing the calculations of various authors 
using various values of b to calculate the cross section. 

2. SELECTI ON OF PHOTONUCLEAR REACTION CROSS SECTION DATA 

Following references 1 and 2, o(y, jn) a is defined as the cross section in which j 
and only i neutrons are emitted and 

o(y, n) = c(y, ln) a ( 2 ) 

by definition. 

Also 

a(y, In) = a(y, n) + a(y, pn) + ••• (3) 

a(y, 2n) = a(y, 2n) a + c(y, p2n) + - (4) 

where the + ••• indicate unimportant additional contributions (for present purposes). The 
total photoneutron cross section is 

c(y, n t ) = o(y. In) + c(y, 2n) + ••• (5) 
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In nucleus-nucleus collisions what is typically measured is 197 Au (RHI, X) 196 Au for 
example, implying that a(y, n) only is needed. Unfortunately what photonuclear 
experimentalist usually measure is G(y, In) or o(y, n t ). Thus we must discuss how to 
arrive at o(y, n) alone. Furthermore, many photonuclear experiments provide contradictory 
data. For this reason I have relied heavily on references 26-28 and have followed their 
recommendations in the selection of data which is summarized in Table 2. 

How to obtain individual (y, n) and (y, p) data for various nuclei is now discussed 
19 7 Au(y. n) 

Following the suggestion of Berman et al 27 ) I have used the data of reference 29 
but multiplied it by a factor of 0.93 (the data was actually extrapolated out to where the 
cross section is zero). The data actually used (Fig. 2, ref. 29) represents a(y. In), but 
because of the large Coulomb barrier, c(y, In) will equal a(y, n). 

89 Y(%-0) 

The data of references 30 and 31 for 89 Y are somewhat different requiring separate 
calculations for each set of data. The data is taken from Fig. 8b of reference 30 and Fig. 3 
of reference 31, both of which represent a(y, In). As for 197 Au, a(y, In) is equal to 
o(y, n) for 89 Y. Following the suggestion of Berman et al 27 ), the data of reference 3 1 has 
been multiplied by 0.82. 

59 Co( y. n ) 

The data of Fig. 3b, reference 32 is used and again o(y. In) is measured but is very 
nearly equal a(y, n) for 59 Co. 

16 Q(y, n ) 

The normalized data presented by Fuller 28 ) represents c(y, ntot)- However this is 
very close to a(y. In) because o(y, 2n) is only about 2% of o(y. In) beyond 30 MeV. 

(This can be seen from Fig. 19A (d) of reference 26). The data of Fuller 28 ) extends out to 
37 MeV and, in principle there is data beyond this that should be estimated and included. 
However various measurements of this higher energy data are not in agreement with the 
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normalized data of Fuller, so it was decided to simply only use Fuller's data up to 37 MeV. 
The neglect of the higher energy photonuclear cross section is compensated by the fact that 
a(y, In) includes not only o(y, n) but also a small component of a(y, pn) which we have 
not subtracted. Further, N(E) is quite small in this energy region and so the error in the 
resulting EM nucleus-nucleus cross section will not be more than a few percent. 

18 Q(y. n't 

Here the data of reference 33, Fig 3b, is used which is for o( y, In). This cross 
section is used for o(y, n) in the present work with no correction. However, based on the 
above discussion for 16 0, the (y, np) correction is expected to be very small. Furthermore, 
the nucleus-nucleus calculation for 18 0 presented herein is only done at 1.7 GeV/N where 
the photon spectrum dominates at low energy, so that a small uncertainty at higher energy, 
in the region of a(y, np) will again not affect the results. 

1 2 C(y. n) 

The data of Fig. 2. 1 reference 28 is used here which represents a(y, ntot)- This 
normalized data represents the original data of Fu ltz et al 34) multiplied by 1.17 which has 
been used to get the data from 30 MeV to 37.5 MeV. Dietrich and Berman 26 ) point out 
that a( y, 2n) is measured to be consistent with zero so that o(y, n to t) ~ cj(y, In). Again I 
have no way of subtracting o(y, np) but this will not affect the results presented herein for 
the same reason as discussed for 18 0(y, n) and 1(i O(y, n) above. 

16 Q(y»p) 

I have used a(y, pt) from Fig. 4.4 reference 28 up to 30 MeV. The data missing 
beyond 30 MeV affects the results of the present work by only a few percent. Following 
the discussion above comparing a(y, 2n) to a(y, nj) for 16 0, the contribution of a(y, 2p) to 
a(y, pt) is expected to be negligible, as is the contribution of a(y, np). 

18 0 (£ P) 

This is given directly in Fig. 3a of reference 33. 
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12 C(y, p) 

Fig. 2.1 of reference 28 gives a(y, Pt) up to 30 MeV. The same considerations for 
16 0(y, p) were followed for 12 C(y, p). 
iS Qfv. M 

This is given directly in Fig. 3c of reference 33. 

197 Au (y. 2n) 

As for 197 Au (y, n) (see above) the data of reference 29, multiplied by a factor of 0.93, 

was used. 

59 Co (y. 2n) 

This data was taken from reference 32. 

Table 2 Cho ice of Photonuclear Reaction Cross Sections 

(see text for details) 


Nucleus 

Data 

Remarks 

* 97 Au 

ref. 27, 29 

Data of ref. 29 are multiplied by 0.93 following 
suggestions of ref. 27. 

89y 

ref. 30, 31 

Data for these two references differ. Thus two sets of 
calculations for both data sets is performed. Data of 
ref. 31 are multiplied by 0.82 following suggestions of 
ref. 27. 

59Co 

ref. 32 

The only existing data. 

18q 

ref. 33 

The most accurate data that exists. 

160 

ref. 28 

The normalized data presented in ref. 28 is used. 

12C 

ref. 28, 34 

The normalized data presented in ref. 28 is used. 
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3. CRITERIA FOR EVALUATING THEORETICAL ANALYSES 


The only adjustable parameter that appears in the present theory is the minimum 
impact parameter b, below which the reaction proceeds via the Strong Interaction. This has 
been the subject of much discussion and every author chooses their own form. For 
instance Heckman and Lindstrom 15 ) and Olson et al 13 ) choose a form 

b = Ro i (p) + Ro l (T) -d (6) 

where Ro i is the 10 percent charge density radii of the projectile and target and d is a 
parameter measuring the amount of overlap. These authors choose values of d ranging 
from 0 up to 3 fm. Hill et al 14 > 20 - 23 ) choose 

b = r 0 [Apl/3 + A t 1/3 - XfAp' 1 / 3 + A T ' 1/3 )] (7) 

where X = 0.75 and r 0 = 1.34 fm. Further, Bertulani and Baur 25) use the form 

b = R(p) + R(T) + nil a (8) 

Z 1 Z 2 e " 

where a is given by — with m 0 the reduced mass and v the relative speed. The 

m v 
o 

above authors variously suggest that a particular form of b accounts for Rutherford bending 
of the orbit (derivation from a straight line) and the effect of a finite charge distribution. 

The problem with choosing a particular form of b and then comparing a resultant 
theoretical cross section to experiment is that the theory incorporates (perhaps unjustified) 
assumptions concerning b. Then, when one compares theory to experiment and makes 
claims about the WW method's validity, it is not only the WW method that one is testing 
but also mixed in is a test of one’s assumption for b. I believe that this approach which has 
been taken in the literature leads to ambiguity concerning whether the WW method agrees 
with experiment. 

Let us now formulate a model-independent criterion that will enable us to firmly 
establish whether or not WW theory agrees with experiment. The simplest possible 
assumption that one can make concerning b is that it is the sum of the projectile and target 
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charge radii equivalent to choosing d = 0 in equation 6 or X = 0 in equation 17 or a = 0 in 
equation 8. All cross section calculations listed in Table 3 are calculated using this naive 
assumption for b. If theory agrees within experimental error, then this is taken as 
indicating that WW theory agrees with experiment. However, what is also calculated is the 
value of d from equation 6 needed to make theory agree with experiment. In this case d 
simply determines the difference of b from our naive assumption of the sum of the radii. 
Remember that the theoretical cross section o is calculated for d = 0. If agreement between 
theory and experiment is found it simply mean's that the value of bmin (or d) needed to fit 
experiment is the sum of the radii (or d = 0). In this case d has a reasonable value and it 
can be claimed that WW theory agrees with experiment. If the theoretical cross section (for 
d = 0) does not fit experiment, but does fit it for another reasonable value of d (say 0.5 fm) 
then again agreement with WW theory and experiment is claimed given the inherent 
uncertainties in b. On the other hand if the calculated cross section (for d = 0) does not 
agree with experiment and if a ridiculously large value of d (say 10 fm) is needed for 
agreement then we conclude that WW theory does nsl agree with experiment. Finally if an 
intermediate value of d (say 3 fm) is needed for agreement then the validity of WW theory 
is uncertain. These criteria however are best expressed in terms of percentages. The 
percentages listed in Table 3 for d are the percentage of d relative to Ro i(p) + Ro-l(T). 
Thus the criteria are re-stated as follows: 

a) Where a zero value of d is listed, then WW theory does agree with experiment. 

b) Where d is say 30% or greater of R 0 -i (p) + R o l CO then WW theory 
definitely does not agree with experiment. 

c) Where d is between 0% and 30%, then a definitive statement concerning 
agreement or disagreement is not possible. 

Note that this 30% criteria was chosen to be as pessimistic as possible. One could argue 
that it should in fact be lowered. However, its main use is to highlight any disagreement 
between theory and experiment. This 30% criterion is based on the assumption that 
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Rutherford bending or finite charge effects can not account for differences in d larger than 
30%. 

The advantage of the above criteria is that we can make a model independent 
definitive statement concerning whether WW theory is adequate. Further, we can also say 
where it is clearly inadequate. This latter point means that a clear delineation is possible of 
where any new physics may emerge. The adoption of the above criteria is suggested in all 
future analysis both for WW theory and other theories of N(E). 

4. RESULTS AND DISCUSSION 

The basic calculational method has been described and criteria established for evaluating 
theoretical comparisons. Ten percent charge radii are determined from the compilations of 
references 35 and 36 and are listed in Table 3. For the sake of comparison the theoretical cross 
section has been calculated using the sum of the projectile and target radii as the minimum impact 
parameter. However, as emphasized above, one should not focus on the calculated versus 
experimental cross section, but rather on the minimum impact parameter b needed to fit the data 
which is listed in the second-to-last column in Table 3. Also listed is the difference between the 
fitted b and the sum of the nuclear radii. This difference is denoted by d (last 2 columns Table 
3) as given in equation 6. Note that the values of the fitted b and d are not unique because of the 
non-zero size of the experimental errors. However, where d is listed as zero, it means rather that 
the calculated o based on b as the sum of the radii already agrees with the experimental a (In this 
case no fitting of b or d took place), d is also listed as a percentage of Ro i (P) + Ro-1 (T). 

Where a zero value of d is listed, I conclude that WW theory agrees with the data 
within the experimental uncertainty. It can be seen that this is the case with all the data on 
12 C and 160 projectile breakup. 15 ) This is in agreement with the original conclusions of 
Heckman and Lindstrom I 5 ), although more accurate data could conceivably change this 
situation. With the data for 18 0 projectile breakup and 197 Au, 89 Y and 59 Co target breakup 
the situation is more complicated. First of all, note that even though the experimental and 


12 



calculated cross section might only differ by a small amount, the value of the overlap 
parameter d required to fit the data may be enormous. This is another reason why one 
should concentrate on d and not a in comparing WW theory to experiment. 

Based on the above 30% criterion, I conclude from Table 3 that WW theory 
disagrees with experiment for the following reactions: 


+ Ti !?N + X 

(1.7 GeV/N) 

12c + !97 Au -» 196 Au + X 

(2.1 GeV/N) 

20Ne + 197 Au -> !96Au + X 

(2.1 GeV/N) 

20 Ne + i97 Au — » 195 Au + X 

(2.1 GeV/N) 

40 Ar + J97 Au — » 195 Au + X 

(1.8 GeV/N) 

20Ne + 89y -» 88 Y + X 

(2.1 GeV/N) 

40 Ar + 89y 88y + X 

(1.8 GeV/N) 

12C + 59 Co -> 57 Co + X 

(2.1 GeV/N) 

16o + i97 Au -> i96 Au + X 

(60 GeV/N) 

16q + l 97 Au l % Au + X 

(200 GeV/N) 


If one lowers the criterion to say 25% then WW theory disagrees with the following 
additional reactions: 

180 + Ti-»17 0 + X (1.7 GeV/N) 

l8o + 23«U -* 17 0 + X (1.7 GeV/N) 

It can be seen that there is serious disagreement between WW theory and 
experiment for l 97 Au target fragmentation both at low (2.1 GeV/N) and high energies (60 
and 200 GeV/N), for both single and double neutron emission. There is also disagreement 
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for two of the 89 Y data sets. (The disagreement for 18 0 — » 17 N and 59 Co — » 57 Co is 
perhaps due to the fact that the measured and calculated cross sections are very small.) It is 
also likely (using 25% criterion) that WW theory does not account for some of the 18 0 
projectile fragmentation data at 1.7 GeV/N. 

Summary and Conclusions 

By considering the value of the overlap parameter d, needed to fit the data I have 
been able to clearly delineate the region where WW theory disagrees with experiment. My 
basic conclusion is that target fragmentation of 197 Au is not understood either at high or 
low energy. This is very important. The results of Hill et al 20-23) did indicate that WW 
theory failed for 197 Au at 60 and 200 GeV/N and thus one would naturally conclude that 
the failure is due to a high energy effect. (These authors plotted the cross section as a 
function of nuclear charge to see if it followed the trend of WW theory.) However, using 
the present criterion based on the 30% overlap parameter it is clear that the failure for 197 Au 
also occurs at low energy (2.1 GeV/N). Thus I conclude that the problem may be with the 
nature of 197 Au EM fragmentation and not necessarily due to high energy assumptions. 
(Note that the photonuclear data for 197 Au is extremely accurate. 27,29)) High energy data 
for lighter nuclei (say 59 Co) would refute or verify this conclusion. 

A referee has pointed out that for heavy nuclei such as 197 Au and 238u, the first 
order perturbation theory on which the WW theory is based is probably incorrect and that 
the WW formula is only the first approximation. This could explain why WW theory fails 
for 197 Au and suggests that it would be worthwhile to calculate higher order effects. Note 
that 197 Au is the heaviest nucleus considered in the present work. 

Finally, the whole situation with respect to EM excitations could be very much 
clarified if we had a clear and precise way of calculating the overlap parameter d. Then one 
could use the actual cross section to compare theory and experiment. 
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Table 3 Overlap parameters d required to fit data. Zero overlap corresponds to a minimum impact parameter equal to the sum of 

the nuclear radii. Where d is listed as 0, it means that c calculated for b = Ro i (p) + Rq.i (T), already agrees within error. Data 
are from references 13, 15, 20, 21, 22, 23. 
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